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FOREWARD 


This document is a supplement to the Design Notebook for 
Naval Air Defense Simulation* SSP-G00-06-0049-82 dated 15 September 
1982fVThe supplement Isa detailed description of the AAW Surveillance 
Simulation which was originally developed for OP-654 and was expanded 
under the current contract to model surveillance by radar as well as 
ELINT satellites of naval and air forces. Base watching, various forms 
of barriers and flight following tactics have been provided. Space 
based surveillance vehicles are limited to circular orbits._ 
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INTRODUCTION 


TRWs Waterwheel Program office has been developing the COVERAGE 
simulation for OP-654, OP-96, other customers, and under company IR&D funding 
for several years. Its current purpose and capability Is to simulate 
detections of aircraft, by active and passive satellites and OTH radars. The 
COVERAGE simulation requires several additional programs as pre-processors and 
post-processors. The purpose of this documentation Is to describe Its 
capabilities and purposes. Its relationship to other programs, and to give a 
sufficiently detailed description of the code to permit computer programmers 
to modify the code. 

At this writing, the COVERAGE simulation Is undergoing modifications 
as a part of current Waterwheel studies. In some cases, especially for 
subroutines BRBLD and BRRIER, the documentation describes modifications about 
to be made to the code rather than the current code which will shortly be 
obsolete. 
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OVERVIEW 
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1.1 BRIEF DESCRIPTION OF COVERAGE SIMULATION 

(U) COVERAGE Is a time-stepped Monte Carlo simulation which models 
detections of aircraft and ships by active and passive satellites. A high 
level flow diagram Is shown In Figure 1-1. Satellite and target (i.e., 
aircraft and ship) motion Is modeled explicitly. At each time step the 
position of the satellites and targets are updated and tests are performed 
to see If any of the satellites detect any of the targets. The tests for 
passive detection are: 

o is the target emitting? 

o does the emitter Illuminate any spacecraft? 

o is a ground station visible to the Illuminated spacecraft? 

The tests for active detection are: 

o is the target emitting? 

o does the emitter Illuminate any spacecraft? 

o Is a ground station visible to the illuminated spacecraft? 

o Is the target In the spacecraft's detection annulus (of which 
the outer edge Is power limited and the inner edge is clutter 
limited)? 

o for aircraft targets, does the doppler return from the aircraft 
exceed the return from the background clutter? 

o Is the radar beam pointed at the target? (Note that this 
simulation models mechanical scanning only and not phased 
array.) 

o a random draw against the probability of detection. 

Active and passive satellites can be modeled simultaneously so that, for 
example, an active satellite can be cued to look at a target that was 
Initially detected by a passive satellite. 

(U) In each replication of a Monte Carlo run, the targets commence 
their courses at a random time, so that the satellites are In a different 
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Figure 1-1. (U) Flowchart for COVERAGE (U) (Page 1 of 2) 





























































position when the bombers take off and when ships arrive In an area of 
Interest. 

(U) The simulation requires three Input files and generates three 
types of output (see Sections 1.3 and 1.4). The Input files describe: 

o satellite and target characteristics (COVIN) 

o target courses (CRS) 

o density of background aircraft traffic (BGTFIL) 

The output files describe: 

o time of detection and target position when detected (OBS, 

DETIM, SCAN, COVECH) 

o satellite hardware performance (SCAN and COVHST) 

o data for input to a graphics program (SCPLT, TGPLT, OVPLT). 

(U) COVERAGE is used In conjunction with other programs (see 
Figure 1-2). They are: 

PROFIL - generates CRS Input files for COVERAGE and AAWSIM 

CRSGEN - creates a graphics input file from a CRS file 

STATCV - processes DETIM to generate statistics on coverage 

gaps 

WRNCOV - processes OBS to compute probability of Initial 

detection as a function of target distance along 
course 

AAWSIM - a one-on-one Monte Carlo simulation using CRS and 

OBS to model engagements between an Interceptor 
assisted by early warning aircraft and a bomber in 
which the Interceptor and early warning aircraft 
are launched and vectored based on external 
surveillance reports. 

MAPPRO - a graphics program using DISSPLA subroutines used 

to draw snapshots of COVERAGE and AAWSIM 
scenarios. 

(U) All these programs are written In FORTRAN. The most current 
versions of all these programs compile, load, and execute on CDC computers 
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Relationship of Coverage to Other Programs (U) 




























and on DEC VAX/780 VMS computers. Earlier versions of COVERAGE, AAWSIM, 
and COVSTAT have compiled, loaded, and executed on a MODCOMP Classic III 
computer. 

(U) These programs have been under development for several years 
with funding from: 

o the Systems Analysis Division of the Chief of Naval Operations 
Program Planning Office (OP-96) 

o the Strategic and Theatre Nuclear Warfare Division of the Chief 
of Naval Operations Plans, Policy and Operations Office (OP-65) 
o contractor IR&D funds 

1.2 PURPOSE OF THE COVERAGE SIMULATION 

(U) The overall purpose of the COVERAGE simulation and related 
programs is to assess the utility of candidate surveillance architectures 
to Naval and Air Force anti-air warfare, primarily against the Soviet 
Backfire bomber, and to naval anti-surface warfare, primarily to support 
over-the-horizon targeting of the Tomahawk anti-ship missile. Production 
runs can be made to study mission utility as a function of: 

o the number of satellites 

o satellite orbital parameters 

o satellite performance characteristics 

o number and location of ground stations 

o C 3 time late 

o tradeoffs between satellites and non-space surveillance systems 

o barrier placement (active satellites only) 

o tasking (active satellites only) 

(U) Passive satellite performance Is modeled In terms of the 
spacecraft-target geometry required for detection. Active satellite 
performance is modeled in terms of slew angles and antenna feed motion 
required for mechanical scanning, energy balance, and the number of 
targets that can be tracked simultaneously. 

(U) Non-space surveillance systems that can be modeled (and the 
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program that models them) are: 


o Early warning aircraft, e.g., the Navy E-2C and the Air Force 
E-3A (AAWSIM) 

o Interceptors, e.g., F-14 and F-15 (AAWSIM) 

o High frequency Over-the-Horizon Backscatter (OTH B) radar (to 
be added to COVERAGE) 

o DEW Line radars (to be added to COVERAGE). 

(U) Tasking of active satellites addresses how many carrier battle 
groups and how many bomber approaches can be adequately covered 
simultaneously. 

(U) Mission utility for anti-air warfare, or more specifically. 
Naval Anti-air Warfare (AAW) and Air Force Airborne Tactical Warning (ATW) 
to defend North America, is quantified in terms of 

o the probability of initial detection of a raid as a function of 
distance from the carrier battle group or from CONUS targets 
(quantified using COVWRN) 

o the probability that an Interceptor can visually identify 

and/or shoot at a bomber beyond some keepout range (quantified 
using AAWSIM) 

o coverage gap statistics (quantified using COVSTAT) 

Mission Utility for Naval Anti-Surface Warfare can be quantified, 
currently, in terms of coverage gap statistics and probability of Initial 
detection. 
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1.3 


INPUT 


1.3.1 BGT FIL Input File 

The BGT FIL Is a file describing worldwide background aircraft traffic 
densities. There are 24 records, each record giving the densities at 000 HRS 
GMT, 0100 HRS GMT, ....,2300 HRS GMT for 10° x 10° latitude-longitude bins. 
Densities are In number of aircraft per million square miles. Each record has 
17 rows corresponding to 10° latitude increments from ± 90° to - 90° and 36 
columns corresponding to 10* longitude bins from 0° to 360°• 

OBS Output File 

The observation output file is written on (unit 12) every time a 
sensor detects an aircraft target. Data written are: 

Column 1 Time of detection (hours) since simulation start; 

Column 2 Number of sensor making the detection; 

Column 3 Course number; 

Column 4 Attack (replication) number; 

Column 5 Target number (row on TARGTS array); 

Column 6 Target latitude (degrees); 

Column 7 Target longitude (degrees); 

Column 8 Time (hours) since target began course; 

Column 9 Bearing in degrees (relative to north) from target to 

sensor; 

Column 10 Distance (n.ml.) of target from destination; 

Column 11 1 - target Is classified; 0 - target Is not 

classified; 






1.3.2 COVIN Input File 


NATTAK 

Integer 

- 

Number of replications; 

NCRS 

Integer 


Number of aircraft and ship 
courses to be read in from CRS 
file; 

LSEED 

Integer 

Odd number 

Random number seed; 

BRTST 

Logical 

Always FALSE 

Obsolete; 

SIMARR 

Logical 

TRUE 

FALSE 

All bombers arrive at missile 
launch points simultaneously; 

All bombers take off 
simultaneously; 

TLIH1 

Real 

Hours 

(Passive sensor only) Emission 
duration immediately after 
aircrafts take off; 

TLIM2 

Real 

Hours 

(Passive sensor only) Aircraft 
emission duration Immediately 
prior to missile launch; 

Rl, R2 

Real 

Hours 

(Passive sensor only) Aircraft 
emitter usage ratio between 
emissions at beginning and end of 
flight. Emitters are used Rl 
hours out of every R2 hours; 

THRSH 

Real 

Knots 

(Active sensor only) Minimum 
discernable velocity; 

BELEV 

Real 

Degrees 

(Active sensor only) Outer 
annulus edge at which aircraft 
barriers are searched (earth 
grazing angle); 

BW 

Real 

Degrees 

(Active sensor only) Beamwldth; 

PD 

Real 

0-10 

(Active sensor only) Probability 
of detection for single dwell; 

PC 

Real 

0 - 10 

(Active sensor only) Probability 
of classification - obsolete; 

TPRINT 

Real 

Hours 

Obsolete; 

TPRINT2 

Real 

Hours 

Obsolete: 





IFLG1 

Integer 

0 

Do not compute rotation of node; 



¥ 0 

Compute rotation of node; 

I HE AO 

Integer 

0 

Suppresses Input echo on unit 6; 



1 

Prints input echo on unit 6; 

I PLOT 

Integer 

0 

Suppresses output on SCPLT and 
TGPLT files (in main program); 



1 

Prints output on SCPLT and TGPLT 
files (in main program); 

IPLOT2 

Integer 

0 

Suppresses output on DETIM file 
(in SETL) 



1 

Prints output on DETIM file (in 
SETL); 

LAA 

Logical 

TRUE 

(Active sensor only) Debug output 
on SRT, AA, and CLM arrays; 



FALSE 

Suppresses debug output; 

NSCAN 

Integer 


Replication for which SCAN, CVPLT, 
SCPLT, and TGPLT output Is 
desired. Zero If no output Is 
desired for these files; 

TPRINT 

Real 

- 

Obsolete 

NSW 

Integer 

1 

Satellite coverage area Is a 
circle (one edge); 



2 

Satellite coverage area Is an 
annulus (two edges); 

ISW 

Integer 

1 

Method of Inputlng satellite 
coverage area edges; 

Earth central angle for each edge; 



2 

Earth grazing angle for each edge; 



3 

Spacecraft look angle (from nadir) 
for each edge; 

TMAX 

Real 

Hours 

Time for each replication or 
attack; should be longer than the 
longest course time; 

NSC1 

Integer 

1-8 

Number of type 1 sensors; 

LTYP1 

Logical 

TRUE 

Type 1 sensors are active; 



FALSE 

Type 1 sensors are passive; 








DTIM1 Real Hours 

OESC 


0ESC(1, J) 

Real 

Feet 

0ESC(2, J) 

Real 

0.-1.0 

0ESC(3, J) 

Real 

Degrees 

0ESC(4, J) 

Real 

Degrees 

OESC(5, J) 

Real 

Degrees 

0ESC(6, J) 

Real 

Degrees 

0ESC(1, J) 

Real 

Degrees 

0ESC(2, J) 

Real 

Degrees 

OESC(3, J) 

Real 

Oegrees 


OESC(4, J) 

Real 

Degrees 

OESC(5, J) 

Real 

n.ml. 

0ESC(6, J) 

Real 

n.ml. 

ORBPAR(NSC, 
1-2, ISW) 

Real 

Degrees 

BATCAP 

Real 

Kilowatt Hrs 
(KWH) 

PSOLAR 

Real 

KW 

PSPC 

Real 

KW 


PPL ON 

Real 

KW 

PPLOFF 

Real 

KW 


Time Increment for Type 1 
spacecraft; 

Orbital parameters for spacecraft 
or dimensions of coverage area for 
OTH radars and DEW Line; 

Orbital parameters (initial 
conditions); 

Length of semi-minor axis; 
Eccentricity; 

Inclination; 

Longitude of ascending node; 
Argument of perigee; 

Mean anomaly; 

OTH B/DEW Line Dimensions; 

Latitude of coverage center; 
Longitude of coverage center; 
Heading (from north) from center 
of coverage area axis of symmetry 
(or bisector); 

Half angle of width of coverage 
area; 

Radius to Inner edge; 

Radius to outer edge; 

Satellite coverage area dimensions 
(See ISW); 

(active satellite only) Maximum 
storage capacity of battery; 
(active satellite only) Average 
n power Input due to Insolation; 
(active satellite only) Average 
power required for base load of 
spacecraft; 

(active satellite only) Power 
required with payload on; 

(Active satellite only) Power 
requled with payload off; 





CLSWID 

Real 

Degrees 

(active sensor only) Maximum feed 
angle from boresight. (rotational 
degrees about suborbital point); 

TGTTRN 

Real 

Degrees 

Maximum expected aircraft turn; 

TGTSPD 

Real 

Knots 

Maximum expected target speed 
change; 

DCPKPW 



Obsolete; 

DUTCYL 



Obsolete; 

SRCHRT(3) 

Real 

Degrees/Minute 

(active sensor only) Search rate 



for barrier search for: 

1 Aircraft barrier at outer edge 
of annulus; 

2 Aircraft barrier at inner 
edge of annulus; 

3 Ship barrier at outer edge 
only; 


TTURN 

Real 

Minutes 

(active sensor only) Feed turn 
time during search; 

TRKCYL 

Real 

Hours 

(active sensor only) Track cycle 
time (minimizes energy usage); 

EFF 



Obsolete; 

REQPEN 

Real 

Ratio 

Feed turn penalty for 
reacquisition; 

TIMTRK 



Obsolete; 

PURTIM 



Obsolete; 

NSCZ 

Integer 

2-20 

Number of type 2 spacecraft; 

LTYPZ 

Logical 

FALSE 

Type 2 spacecraft must be passive; 

DTIMZ 

Real 

Hrs 

Time step for type 2 spacecraft; 

NSMCON 

Integer 

- 

Number of ground station sets; 

ISMTAB 

Ground station set specifications; 

ISMTAB(1, J) 

Integer 

1-15 

Set Number; 

ISMTAB(2, J) 

Integer 

- 

Target one start number; 

ISMTA8(3, J) 

Integer 

- 

Target one end number; 

ISMTAB(4, J) 

Integer 

- 

Target two start number 


(first 6/S); 
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ISMTAB(5, J) Integer 


Target two end number 
(last G/S); 

ISMTAB(6, 0) Integer - Spacecraft start number (0 = any 

S/C); 

ISMTAB(7, J) Integer - Spacecraft end number (0 - any 

S/C); 


Notes on ISMTAB; 

1. Target numbers refer to rows on TARGTS array; 

2. Ground stations should be entered sequent!lally at the end 
of the TARGTS array; 

3. Entries 1-3 are used for passive sensors only and; 

4. If active satellites are modelled, ISMTAB refers to active 
satellite ground stations only; 


ITSTOR 

Integer 

Zero 

Never use; 

NUMCT 

Integer 

Always 1 

Number of simultaneous observation 

NTRGTS 

Integer 

1-35 

by spacecraft required for 
detection; 

Number of targets to be Input to 

TARGTS 

Real 

the TARGTS array; 

Input area for; 

o Aircraft 

o Ships 

o Aircraft barrier; 
o Ship barrier; 
o Ground stations; 


TARGTS 
TARGT$(J, 1) 

TARGTS(J, ?.) 
TARGTS(J, 3) 
TARGTS(J, 4) 
TARGTS(J, 5) 


Input for aircraft and ships; 

1 Aircraft raid count; 

1 Ships raid count; 

Course number; 

(Passive sensor only) Emitter scan width In azimuth; 

(Passive sensor only) Emitter elevation limits 1 or 2; 

(Passive sensor only) Emitter's lower elevation limit (grazing 
angle degrees above local horizon); 




TARGTS(J, 6) 

TARGTS(J, 7) 
TARGTS(J, 8) 

TARGT$(J, 9) 


(Passive sensor only) Emitter's upper elevation limit (grazing 
angle degrees above local horizon); 

Course start time (hours) - computed by program. Input a zero; 
For passive sensors. Input time (hours) after takeoff at which 
emitter first turns on; 

For active sensors, enter any negative number; it will be 
replaced with the time of initial detection; 

Time (in hours) of course duration; is computed by program; 
input a zero; 


Aircraft and Ship Barriers; 

Note; Ship barriers must be circles centered at a carrier battle group; 
aircraft barriers can be circles or annular sections centered anywhere; 


TARGTS(J, 1) 
TARGTS(J, 2) 
TARGTS(J, 3) 
TARGTS(J, 4) 

TARGTS(J, 5) 
TARGTS(J, 6) 

TARGTS(0, 7) 


TARGTS(J, 8) 
TARGTS(J, 9) 


Must be zero; 

Latitude (degrees) of barrier center; 

Longitude (degrees) of barrier center; 

Inner radius (n.mi.) of barrier (must be zero for 
circular barrier); 

Outer radius (n.mi.) of barrier; 

Bearing (degrees from north) of annular barrier 
axis of symmetry; 

Barrier ID. Must be an odd number from 1 to 19 for 
aircraft barrier and any integer from 21 to 30 for 
ship barriers; 

Barrier width (degrees); 

Barrier type (1 for aircraft, 2 for ships); 


Ground Stations: 


TARGTS(J, 1) 
TARGTS(J, 2) 
TARGTS(J, 3) 
TARGTS(J, 4) 
TARGTS(0, 5) 


Must be zero; 

Latitude (degrees); 

Longitude (degrees); 

Elevation coverage limits (1 or 2); 

Earth Grazing angle (degrees) of lower elevation 
limit); 
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TARGTS(J, 6) 


Earth Grazing angle (degrees) of upper elevation 
limit; 

TARGTS(J, 7) ID; Must be any negative Integer (entered as a 

real); 

TARGT$(J, 8-9) Not used; 

Revision data: 


I 

Integer 

< NSC1 

Number of type 1 sensors to be 
modelled; 

J 

Integer 

£ NSC 2 

Number of type 2 spacecraft to be 
modelled; 

IDO 

Integer 

£ NSMCON 

Number of ground station sets to 
be modelled; 

K 

Integer 

< NTARGTS + 

NSMCON 

Number targets to be modelled; 





1.3.3 COURSE FILE 

(U.) This file defines all bomber flight paths to be simulated in terms 
of state vectors. Each state vector consists of: 

1) Time of course change (hours) 

2) Lattitude of the course change (degrees) 

3) Longitude of the course change (degrees) 

4) Target acceleration after course change (knots/hr.) 

5) Target velocity after course change (knots) 

6) Time since last course change (hours) 

7) Distance flown since last course change 
(nautical miles) 

8) Total distance flown thus far (nautical miles) 

9) Target heading after course change (degrees) 

, 10) Distance from carrier (nautical miles) 

A separate program (PROFILE) generates elements one through nine for each the 
state vectors. AAWSIM computes the tenth element of each state vector. 
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1.4 


OUTPUT 


1.4.1 DETIM Output File 

OETIM ts written on (unit 8) by subroutine SETL every time a new 
contact (or detection) begins of a barrier, aircraft, or (for passive sensors 
only) a ground station. Data written are: 

Column 1 Target number (rows or TARGTS array); 

Column 2 Attack/replication number; 

Column 3 Time (hours) since aircraft began course or since (for 

ground stations and barriers) replication began; 

Column 4 Down or gap time (hours) since end of last contact and 

start of this contact; 

Column 5 0 - target has not been classified; 

1 - target is classified; 











1.4.2 TGPLT, SCPLT, and CVPLT Output Files 


These files have a special format required by the MAPPRO graphics 
program. This format is 


FILID 
data rows 
0,0,0,0 
FILID 

data rows repeated an unlimited number of times 

0,0,0,0 

End-of-file 


The format for each data row is LONGITUDE, LATITUDE, symbol number, 
heading where all angles are In degrees. Symbol numbers are the following 
integers: 


1 

0 

- 1 
- 2 

- 3 

- 4 

- 5 

- 6 
- 6 


Causes a DISSPLA symbol to be drawn for each data row; 
A line connects all the points indicated on the data 
rows; 

Do not use; 

Do not use; 

Draws Backfire symbol; 

Draws F-14 symbol; 

Draws CV symbol; 

Draws OTH B symbol; 

Draws other symbols satellites, RPV, E-2C, or an "X"; 


FILIDs must be ascending real numbers. Free format is used. 

TGPLT writes a data row to show the current position of every bomber 
at every simulator time step. The FILID Is the simulator time at which the 
data row Is written. 

SCPLT writes a data row at every time step for 

o suborbital point of each satellite; 
o boreslght almpolnt for each cluster (active 
satellite only) 
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o Intersections of annulus edges with aircraft and 
ship barrier edges; 

Again, the FIL1D Is the simulator time at each time step. 

CVPLT writes scenarlc data rows. Data Is Indicated by the FILID. 

FIL1D DATA 

1-50 Aircraft and ship courses; 

51 - 100 Destination of aircraft and ships; 

101 - 150 Barrier centers; 

200 Aircraft and ship positions when detected; 




1.4.3 SCAN Output Files 


The SCAN file shows the status of each sensor at each time step, 
especially by indicating which targets are observed. Data written for each 


sensor are: 


to? 


Column 1 Time (hours); 

Column 2 For each target In TARGTS array o if not detected; 

9 If aircraft cannot be detected because velocity 
component along line of sight is less than the minimum 
dlscemable velocity (active satellite only); 

7 aircraft (ship) not detected because probability of 
detection is less than one (active satellite only); 

A-Z indicates 

o barrier aircraft, or ship is detectable; 
o cluster ID; 

Column 3 # Indicates a scheduling conflict 

blank Indicates no scheduling conflict (active 
satellite only); 

Column 4 letter Indicates which cluster was scheduled;(The 

remaining data are for active satellites only.) 

Column 5 Boreslght slew angle since last time step (rotation 

angle degrees about boreslght); 

Column 6 Off-boreslght feed angle required to scan cluster 

(rotation angle degrees about suborbital point); 

Column 7 Aircraft track count (the number of aircraft being 

tracked at current time step); 

Column 8 Battery state charge (0 to 100*) ; 


Also, there Is a summary column, which shows for each target. 

8 target Is detected 
0 target Is not observable 

018 target Is observable by I satellites but Is 
not detected because of scheduling conflicts; 
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2.0 OVERVIEW OF MAIN PROGRAM 


This chapter describes the functions performed by the main program 
code. The code Is referred to In terms of the FORTRAN DO/LOOP and other 
labels and subroutines. The main program Is long and performs much of the 
simulation logic. The overall structure of the main program (called COV) 
Is: 


Declaration Statements 
Comments describing Inputs 
Data statements 
FORMAT statements 

Input and Initialization (labels 10-200) 

Time loop (labels 200-1000) 

Sensor loop (labels 295-600) 

Target detection loops (labels 380-500) 


Target detection loops Include logic or subroutine which test for the 
following sensor - target combinations. Under the LOCATION heading, names 
refer to subroutines and numbers refer to FORTRAN labels and DO LOOPs in' 
the main program. 


SENSOR TARGET LOCATION 

active satellites aircraft barriers BARRIER (In LOOP 385) 
aircraft BARRIER, LOOP 500, and 

828-855 

ship barriers BARRIRS 

ships BARRIRS, LOOP 500, and 

828-855) 

ground stations 355-380 

passive satellites aircraft, ships, LOOP 500 

and ground stations 

OTH, B radar & aircraft FXSENS 

DEW/LINE 







Other Important sections of code In COV, and the associated sensor. If 
applicable, are: 


SENSOR 

FUNCTION 

LOCATION 


Calculate start & end 
time for each replication 
and takeoff times for 
aircraft 

200-210 

all satelltes 

Solar eclipsing 

220-230 

ECLPSE 


Update position of ships 
and aircraft at each time 
step 

260-295 

all satellites 

Update spacecraft position 
(In terms of suborbital point) 

300-320 

CORBIT 

active satellites 

Test for ground station 
visibility 

355-380 

active satellites 

Identify candidate search 
areas at which boreslght 
should be pointed 

FSTEP 

500-580, Including 
CLUSTR 

passive satellites 

Test for ground station 
visibility 

LOOP 625 

active satellites 

Decide at which search area 
to point the boreslght of 
each spacecraft 

CLSTPP, just 
after 675 

— 

Output SCAN and plotting files 

900-997 


The following test describe COV sequentially from label 10 to the end In 
greater, but not complete, detail: 


LABEL/LOOP/SUBROUTINE DESCRIPTION 

LOOP 10 Rewind files. 

COMBLK Initialize earth gravitational constant 

(GMU), earth rotation rate (ME) and earth 

radius (RE) In distance, time, angle units of 
feet, hours, and degrees. 


SECOND 


Save current value of CPU time 








ZERRV 

ZERIV 


Set arrays to zero. 


20-30 


40-45 

45-50 


LOOP 75 


LOOP 87 
LOOP 90 


# 


LOOP 100 


105-140 


LOOP 155 


LOOP 165 


165-185 

BRBLD 

BRBLDS 

UPR1 

JD0R6F 

LIHTRP 



Read COVIN file (describes sensors and 
scenarios). 

Echo COVIN file on COVOUT. 

Determine time step (DTIME) and IRATIO, which 
makes It possible toupdate and test the two 
types of satellites at different time steps 
(DTIM1 and DTIM2). TIME Is subsequently set 
to zero. 

Read the CRS file (describes routes followed 
by aircraft and ships) and compute distance 
between starting point of each aircraft or 
ship course and the location of the carrier 
or base that the aircraft or ship Is 
attacking. 

Write aircraft/ship courses on COVPLT file. 

Write locations of carriers or bases that 
aircraft or ships are attacking on COVPLT 
file. 

Echo CRS file on COVSUM file. 

Write header for SCAN file. Contents of 
header Is determined by the number of sensors 
and the number of targets (l.e., aircraft and 
ship barrier, alrcrft, ships, and ground 
stations. 

Spacecraft Initialization. DNODEH Is 
procession of the mode In degrees. XM Is 
spacecraft velocity and HA and HP are apogee 
and perigee. 

Converts some entries of the TARGTS array 
from real to Integer. 

8RBLD and BRBLDS set up arrays that are 
referred to In processing aircraft and 
ship barriers. UPR1 a uniform random 
number generator and Is called with 
JDORGF to compute a random starting date 
for eclipsing. Remaining code Initializes 
variables for computing active satellite 
energy balance. LINTRP Is a linear 
Interpolation routine. 






2-3 





185-200 


200 


200 to 220 


220-230 

BGTINP 

ECLPSE 


235-250 


250-295 


260-263 

BGNCRS 

ENDCRS 


263-275 


275-295 


LOOP 500 


Compute TMAX, the time at which the 

simulation ends. 

Start of time loop, in which sensor and 
target positions are updated and tests are 
made to see if any sensors detect any 
targets. 

For each attack (l.e., replication), compute 
simulated time at which this replication 
begins (SATTAK) and ends (TATTAK) and at what 
time each of the aircraft/ships begins its 
course (TARGTS (NTAR.7)). 

BGTINP updates clock and reads in new 
aircraft background density data If hour 
has changed. ECLPSE checks for solar 
eclipsing of satellite. 

Initialization required each time the time 
loop Is updated. 

Update position of moving targets and DELTIM 
of fixed and moving targets. DELTIM Is the 
difference between the current TIME and the 
time the replication began (for fixed 
targets) or the time the aircraft or ship 

began Its course (for moving targets). 

Check If aircraft/ship has begun or ended 
course. BGNCRS Is called just before an 
aircraft begins Its course so that coverage 
statistics will include the detection or gap 
occurring when the aircraft begins Its 

course. ENDCRS Is called for the same 
purpose just after an aircraft ends Its 
course. Otherwise, If for example, an 

aircraft course took three hours and the 

first detection did not occur until two hours 
Into the flight, that two hour gap would not 
be Inlucded in CVSTAT statistics If BGNCRS 
were not used. 

Interpolation between two state vectors from 
the CRS bit to get aircraft/ship position, 
speed, and heading. 

When target happens within six minutes from a 
state vector, position, course, and speed are 
obtained directly from the state vector 
Instead of from Interpolation. 

Sensor loop. Performs spacecraft, OTH B, and 
DEW Line tests for current time. 
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295-300 

300-320 

CORBIT 

320-350 


355-370 

FXSENS 

LOOP 385 
BRRIER 

LOOP 

LOOP 500 


405-425 


Ships sensor loop when needed if the two 
types of spacecraft are tested at different 
time increments (DTIM1 and DTIM2). 

Update spacecraft suborbital point for 
current time step. 

Compute earth central angle for inner and 
outer edge of spacecraft detection annulus 
based on spacecraft altitude, and satellite 
look angle from nadir (BETA) or earth 
elevation angle above horizon (EPS). 

Test if active satellites have a ground 
station visible. If not, skip rest of sensor 
loop. 

If current sensor Is an OTH B radar or the 
DEWLSC line barrier, test for detection of 
aircraft and then ship to end of sensor loop. 

Test ifany aircraft barriers can be searched 
at current time by current spacecraft 
active. If it Is, store boundaries of search 
area in SRT array for spacecraft scheduling. 

Test if any ship barriers can be searched at 
current time by current spacecraft, if 
ctlve. If so, store search area boundaries 
on SRT array for spacecraft scheduling. 

For active satellites, tests for detection of 
aircraft that have already been initially 
detected, either in a barrier, by a passive 
satellite, by the DEW line, or by an OTH B 
radar. For passive satellites, test for 
detection of a ship or aircraft and for 
ground station visibility. 

Retrieve target position, course, speed and 
DECTIM from TPOS. 

Test Is target Is in annulas by comparing 
earth central angles between suborbital point 
and target and suborbital point and nular 
edges. 


425-440 




445-465 

SETL 


465-475 

JMEQ 

FSTEP 


500 

500-565 


Complete detection tests for passive 
satellites, namely: is aircraft/ship 

emitting? If emitting, is it emitting in the 
direction of the aircraft? 

Emitter on/off usage is described by TUMI, 
TLIM2, Rl, and R2 (see Table ). If tests 
are passed and ground station visibility Is 
not being tested, the detection is written as 
CVSCAN, CVOBS, and SETL is called. SETL 
updates coverage statistics array for current 
replication and may cause write on CVTAP8 
(Input to CVSTAT). If ground station 

visibility is being tested,then target 
observability is Indicated on the IRFTAB 
array. If current target is a ground 
station, then Its observability Is noted In 
ISCAN and IRFTAB and SETL Is called. 

Perform doppler test for active detection 
of aircraft by calling JMEQ. JMEQ computer 
DDOP, the component of aircraft velocity 
along the Hne-of-sight between the target 
and the spacecraft. For detection to occur, 
the absolute value of DDOP must exceed the 
doppler return from clutter, which Is Input 
as THRESH. If the aircraft Is detectable* 
then FSTEP is called to store aircraft 
position for subsequent use in deciding where 
to point the boresight (In subroutines CLUSTR 
and CLSTPP). 

End sensor loop. 

By now, the SRT array holds the spacecraft 
look angle (from radar) and the bearing 
(relative to spacecraft velocity) for each 
target and for each search area end point 
detectable from the current spacecraft. The 
AA array is now computed. It stores the 
suborbital-polnt-centered angle between every 
pair of targets (and search area end points) 
observable by the current spacecraft. This 
Information is used by CLUSTR. If search 
area endpoints exist, then the search area ID 
Is stored on the AA diagonal, which otherwise 
would be zero. 


565-570 


Optional listing of AA for debug purposes. 





570-580 

CLUSTR 


580-630 

SETL 


# 

650-675 


Computes boreslght pointing options for 
current spacecraft by partioning detectable 
targets and search area into clusters. 
Clusters are groups of targets and search 
areas which can be scanned froma fixed 
boresight position by moving the antenna 
feed. After CLUSTR has been called for all 
spacecraft, then CLSTPP is called (after 
label 678) to decide where to aim the 
boresight of each spacecraft. For fixed 
sensors, i.e., OTH B and DEW line, all 
observable targets are assumed to be in one 
cluster and the CLUSTR functions are 
performed in FXSENS. 

Ground station logic for passive spacecraft. 
The ISMTAB array specifies that particular 
spacecraft (ISMTAB(6-7,I)) can report 
detections of particular target 
(ISMTAB(2-3,I)) only if the spacecraft sees a 
target and particular ground station 
locations (ISMTAB(4-5,I)) simultaneously. 
For targets and ground stations, ISMTAB 
refers to rows on the TARGTS/ITRGTS array. 
If detection occurs, SETL is called to update 
coverage statistics and, perhaps, write on 
CVTAP8, JSCAN Is set, and the latitude and 
longitude of the observed target Is written 
on tape 8. 

ITSTOR is never used. 


675-LOOP 815 Decide at which cluster to aim active 

CLUSTR satellite boresight. This process Is called 

scheduling. CLSTPP knows the clusters of all 
spacecraft. 

LOOP 815 Update 'he ISCAN summary (last column) to 

show how many sensors detect each target. 

LOOP 875 Each sensor has now been assigned, or 

scheduled, to look at a particular cluster. 
This loop updates active spacecraft hardware 
requirement parameters (state of change, feed 
angle, lew angle, and number of aircraft to 
be tracked) and determines, based on a random 
draw against an input probability of 
detctlon, whether a particular ship or 
aircraft Is detected. If detection occurs, 
ISCAN Is set, CVOBS Is written on, and SETL 
is called. 
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00875 - 00855 
STATUO 


00855 - 830 
UPRI (TSS) 
RAN (VAX) 


830-855 

SETL 


855-870 

STATUO 


870-875 

ECLPSF 

PWRBAL 

STATUO 


SCSCH (3,IDO) stores the number of the 
scheduled clusters (in terms of a row on the 
CLM/CLM2/CLM3 arrays). STATUD is called 
twice to update histogram statistics. The 
first call updates the number of aircraft 
that must be tracked, which is the sum of the 
background and threat aircraft in search 
areas and In (smaller search) areas where 
threat aircraft are being reacquired and 
tracked. The second call updates the ratio 
of the number of aircraft that must, be 
tracked (track requirements) to the number of 
aircraft that there is time to track (tfack 
capacity). A message is written whenever,!the 
ratio is greater than one. 

For each aircraft or ship target, a uniform 
random draw is made to test against an input 
probability of detection (PD). If the draw 
is greater than PP, detection does not occur 
and ISCAN Is set to -7, which shows up in the 
SCAN printout as a +7. 

If an aircraft is being detected (in a 
search area) for the first time,the time of 
the initial detection is stored in the eighth 
element of the TARGTS/ITRGTS row for that 
target. For detected targets ISCAN is set to 
the number of the scheduled cluster, SETL is 
called, and CVOBS is written on. 

Based on the boresight aim points at the 
current and previous time steps, spacecraft 
slew angle (centered at the suborbital point) 
Is computed and stored In CLBG (5,5, NSC) and 
STATUD Is called to update slew angle 
histogram statistics. STATUO Is also called 
to update histogram statistics on cluster 
width (the suborbital point centered angle 
between the targets or search area corners at 
the edges of a cluster). j 

Update power balance statistics for current 
spcecraft. ECLPSF is called to determine If 
current spacecraft Is eclipsed. In which case 
solar panels are not charged. 


875-900 For targets currently being detected, set 

time of' last detection(TARGTS(NTAR, 11)) to 
current time. 
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90-995 


900-920 

CTOALF 

J 

920-955 

955-985 

985-995 

995-1000 

1000-1050 

STATPT 

SUMMRY 


Write on output files CVSCAN, CVPLT, CVTG, 
and CVSC If DOUT is true. DOUT prevents 
writes at the smaller time step in DTIM1 and 
0TIM2 are different. 

Write on CVSCAN file. CTOALF converts 
positive numbers stored in ISCAN to alphabet 
letters. Different SCAN formats are used 
depending on whether one or two types of 
spacecraft are being simulated. CVSCAN 
shows, for each sensor, which targets are 
detected, the different clusters, which 
pluster Is scheduled, and hardware 
requirements parameters. 

Writes longitude and latitude of current 
position of detected moving targets on CVPLT 
file. 

Writes longitude, latitude, and heading of 
suborbital points of spacecraft on CVSC 
file. If active spacecraft are being 
simulated, the longitude and latitude of 
boresight aim points are written. If active 
spacecraft are simulated and their annulus 
edges Intersect barriers (thus creating 
search areas), the longitudes and latitudes 
of these Intersections, which jre the search 
area endpoints, are also written on SCPLT. 

Write longitude, latitude, and bearing of all 
moving targets at current time step on TGPLT 


End of time loop. If TMAX is exceeded, exit 
loop and end simulation. Otherwise, 
Increment the number of time steps (ITIME), 
time, clock, and date. 

Call final statistics subroutines, write 
endflles, rewind files, and stop. STATPT 
writes histogram dates on tape 4 and SUMMRY 
computes coverage statistics for last 
replication only and writes then on CVECHO. 
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3. SUBROUTINES 

3.1 INPUT 

3.1.1 SUBROUTINE BRBLD 
PURPOSE: 

Reads Input data about aircraft barriers from targets array and stores 
it in BRREF/IBRREF for processing aircraft barriers in BRRIER. Barriers are 
modeled as circles or annular sections. 

INPUTS: 

NTRGTT 

TARGTS (NTAR.l) 

TARGTS (NTAR,2) 

TARGTS (NTAR.3) 

TARGTS (NTAR.4) 

TARGTS (NTAR.5) 

TARGTS (NTAR.6) 

TARGTS (NTAR.7) 

TARGTS (NTAR,8) 

TARGTS (NTAR.9) 


Number of input rows to TARGTS array inputs for 
aircraft barriers; 

Must be zero if row contains barrier information; 
Longitude of barrier center (deg); 

Latitude of barrier center (deg); 

Radius to inner edge of barrier (n.mi.); 

Radius to outer edge of barrier (n.mi.); 

Bearing of bisector of barrier (deg. from north); 
Fixed target ID; for aircraft barriers must be an 
odd number from 1 to 19; 

Barrier width (half angle); 

Not used 


For circular aircraft barriers, TARGTS (NTAR,4) must be zero and the sixth and 
eighth entries are not read. Annular barriers must have two rows on the 
TARGTS array, the first of which must have the inputs described above and the 
second of which must have the first entry set to zero are the seventh set to 
the same barrier ID. 

TARGTS array inputs for aircraft. 

TARGTS (NTAR,1) Must be greater than zero; 

TARGTS (NTAR,2) Course number In CRS input file; 

TARGTS (NTAR.3) Not used 

TARGTS (NTAR,4) Not used 

TARGTS (NTAR,5 } Not used 

TARGTS (NTAR.6) Not used 

TARGTS (NTAR,7) Time after start of replication that bomber 

begins course. Is not an Input but Is computed by 
program; 






TARGTS (NTAR,8) 


Set to any negative number. If and when aircraft 
Is detected In a barrier, this entry Is set to 
the time at which Initial detection occurs. 

TARGTS (NRAT.9) Length of flight (hours); 

OUTPUTS 

NUMBR Number of aircraft barriers in TARGTS array; 

IBRREF (l.IBR) Barrier ID; 

IBRREF (2.IBR) Number of (first) row in TARGTS array with 

Information about barrier center, radii, and 
bearings; 

IBRREF (3,I8R) Number of (second) row In TARGTS array, with 

barrier ID only; 

IBRREF (4,IBR) Latitude of pole of a great circle going through 

the barrier center and defining one side of the 
barrier; 

IBRREF (5,IBR) Longitude of pole first described; 

IBRREF (6,IBR) Latitude of pole of great circle of radial 

defining other side of barrier; 

IBRREF (7, IBR) Longitude of pole just described; 

Has a shorter current detection period than that 
of any of the targets seen by the other cluster; 





METHOD: 


Aircraft barriers are defined as annular sectors: 



This subroutine searches the TARGTS array to build the BRREF array. The poles 
for the radlals are computed by adding 90 to the bearing of a radial and 
finding a point 1/4 of an earth radius away from the center along that 
bearing. If the inner radius Is zero, then a circular barrier is assumed and 
no poles are calculated. The data described here are used In subroutine 
BRRIER. 






3.1.2 SUBROUTINE BRBLDS 


PURPOSE: 

Sets up a reference array. BREFS, for use by subroutine BRRIRS In 
processing ship barriers. Like BRBLD, which is for aircraft barriers, BRBLDS 
gets the data for the reference array by searching the TAROTS array. Ship 
barriers are models as circles centered at a carrier battle group with an 
Input radius specified on the TARGTS array. They are searched only at the 
outer edge of the annulus. Ship barriers are used to model radar satellite 
searches of enemy ships for providing over-the-horizon targeting of Tomahawk. 

INPUTS: 


For ship barriers, TARGTS must have the following Inputs: 


TARGTS (l.NTAR) 
TARGTS (2,NTAR) 
TARGTS (3.NTAR) 
TARGTS (7.NTAR) 


Zero; 

Latitude of barrier center; 

Longitude of barrier center; 

Barrier ID; must be a number from 21 to 30 
(Implies a maximum of ten ship barriers); 


TARGTS (8.NTAR) Radius of barrier (n.ml.); 


To Input ships,' TARGTS must have: 


TARGTS (NTAR.l) 
TARGTS (NTAR.2) 
TARGTS (NTAR.7) 

TARGTS (NTAR.8) 

TARGTS (NTAR.9) 


A negative number; 

Number of course Input In CRS file; 

Time, (In hrs) at which ship begins Its course 
(computed by program); 

Time at which ship Is Initially detected 
(computed by program); 

Duration of ship course (hrs); 


OUTPUTS: 

NUMBRS 

IBREFS (1,IBR) 
IBREFS (2,IBR) 
IBREFS (3,IBR) 
IBREFS (4,IBR) 


Number of ship barriers Input In TARGTS array; 
Barrier ID number; 

Row In TARGTS array describing this barrier; 
Not used 
Not used 
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IBREFS (5,IBR) Not used 

BREFS (6,IBR) Latitude of barrier center (degrees); 

BREFS (7,IBR) Longitude of barrier center (degrees); 

BREFS (8,IBR) Barrier radius (n.ml.)i 

BREFS (9,IBR) Barrier radius (degrees); 

METHOD: 

As part of program Initialization, BRBLDS tests every row of the 
TAR6TS array for a zero In the first entry and a number between 21 and 30 in 
the seventh entry. If these tests succeed, then the row describes a ship 
barrier and the appropriate Information Is put Into the BREFS array. 





3.2 SUBROUTINE CORBIT 

Subroutine CORBIT calculates the suborbital point and the velocity of 
the suborbital point. 

Spacecraft position and velocity in geocentric inertial spherical 
coordinates are translated Into earth geographical coordinates. The elements 
used to describe the circular orbit for an earth satellite are listed in table 
1 . 

CORBIT first calculates the orientation of the orbit plane, that is, 
the longitude of the ascending node at the time of interest, 

ON » TIME(DNODE-WE) + 00 

and the argument of perigee; the rotation of theorbit within the orbital plane 
PM * (XNU * TIME + 0M) MOO 360° 

where MOO Is the remainder function. 

The position of the satellite at the time of Interest thus calculated 
is translated to earth coordinates. 












TABLE 1 


Symbol 

Element 

Definition 

DI 

Inclination 

Angle between orbit 
plane and earth's 
equatorial plane 

XNU 

Mean Angular Motion 

Angular velocity in the 
orbital plane y/a 3 

A 

Distance 

Distance from earth 



center to the satellite 

DNODE 

Precession of 

Angular velocity of the 


ascending node 

orbit plane about the 

line normal to the 
earth's equatorial plane 
and throuqh the origin 
(earth center) 

PM 

Argument of perigee 

The rotation of the 

orbit within the orbital 
plane 

00 

Right ascension of 

Longitude of ascending 


the ascending node 

node at reference time 

OM 

Mean anomaly at epoch 

360° x At/p where P Is 
the period and t Is the 

time difference between 

the epoch and the last 
perigee passage before 
epoch 

To 

Epoch 

Time for which elements 

are specified 

ON 

Longitude of 

Orientation of the 


ascending node 

orbital plane at the 

time of Interest 

WE 

Earth rotation rate 
(degrees/hour) 


GMU 

Gravitational constant 







The coordinate of the point is (PM, 0) in Inertial coordinates 

In earth reference the old North pole is at (0N-90*, 90-DI), the 
coordinates of the point in the new system 0^= 90° 

COS(90 - LAT) = COS(DI) * 0 + SIN(DI)C0S(PM - 90) 

SIN(LAT) * SIN(DI) SIN PM 
LAT = SIN' 1 SIN(01)SIN PM) 

SlN(L0N-90) = SIN(90 - PM) SIN 0°/SIN(90 - LAT) 

COS(LON) * COS(PM)/COS(LAT) 

LON = COS 1 COS (PM )/COS (LAT)| . 















3.3 ACTIVE SENSOR LOGIC 


3.3.1 SUBROUTINE BGT 
PURPOSE: 

This subroutine is responsible for computing the background traffic 
in a specified area. 

INPUT ARGUMENTS: 

AREA, ALAT, ALON, BGTDEN 

Three calling sequence inputs are used to define the area of 
interest, and a fourth variable, the background density array, is provided 
through X8 Common. 

AREA Area of region in square miles; 

ALAT Latitude representative of region in degrees; 

ALON Lonqltude representative of region in degrees; 

BGTDEN Background traffic density data, stored in an 18 x 36; 
array, in aircraft Der million square miles. 

(BGTDEN is contained in COMMON X8); 

OUTPUT ARGUMENTS: TRAFIC 

The aircraft count resultina from the computations of this routine 
are returned via the calling sequence variable TRAFIC. 

TRAFIC Number of aircraft in the region. 


METHOD: 

This routine computes the aircraft count in a region based upon the 
area and the density associated with the region. The region's latitude 
and longitude are converted Into array subscripts by the following 
formulae: 
















* 


LAT = 

LON = 


9 - 1F1X(1.E-6 + ALAT/10.) 

10 - lFIX(l.E-6 + ALAT/10.) 

18 + lFlX(l.E-6 + AL0N/10.) 

19 + 1F1X(1.2-6 + AL0N/10.) 


if ALAT > 0. 
if ALAT < 0. 
if ALON 1 0. 
if ALON < 0. 


The traffic count Is then computed by the following statement, which 
includes the conversion from aircraft per million square miles: 


TRAFIC * BGTDEN(LAT.LON) x AREA/10.E6 
REFERENCES: 

The routine is called from BGTSCN and CLUSTR. 


6 
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3.3.2 Subroutine BGTINP 
PURPOSE: 

This routine is responsible for keeping the background traffic data 
current with the simulation time. The background traffic file B6TFIL 
contains densities by hour for the ten degree square segments covering the 
surface of the globe. As model time progresses, the density data must be 
revised for each hour. 

INPUT ARGUMENTS: 

CLOCK, BGTTIM, BGTDEN 

The routine receives the simulated clock time (24 hour, with 
fractions of an hour) through the calling sequence. It maintains BGTTIM, 
the time associated with the current density data set, and BGTDEN, the 18 
by 36 array of background traffic densities, the density data set. These 
latter two items are input from the BGTFIL on FORTRAN logical unit number 
20 and are maintained in COMMON X8. 

CLOCK 24.00 hour simulation clock time 

BGTTIM 24.00 hour time, in hourly increments, associated with 
BGTDEN 

BGTDEN Background traffic density, stored in an 18 x 36 array, 
in aircraft per million square miles. 

s 1 

OUTPUT ARGUMENTS: ' 

BGTTIM, BGTDEN 

These arguments are tn' current time and the current background 
traffic density data as uDdated from the BGTFIL as necessary. They are 
defined above. 






METHOD: 


The CLOCK time, maintained in the main program, is compared to the 
BGTTIM value. Initially, B6TTIM is set to -10.E20. If CLOCK is greater 
than BGTTIM, the next record for the BGTFIL file (BGTTIM and BGTDEN) is 
read and the times are again compared. If BGTTIM exceeds 24., the file is 
rewound before it is read. 

REFERENCES: 

This routine is called from COV, the main routine, early in the 
time Iood. 







3.3.3 Subroutine BGTSCN 


PURPOSE: 

This routine Is responsible for accumulating the quantity of 
background traffic In a scan area bounded by points A and B along a 
spacecraft annulus of depth DISTDP. 

INPUT ARGUMENTS: 

The following Is a list of Input arguments passed by the argument 

list: 


SCLAT 

SCION 

ALAT 

ALON 

BLAT 

BLON 

S 

LINNER 
01STOP 


Spacecraft's suborbital point latitude In degrees; 
Spacecraft's suborbital point longitude In degrees; 
Latitude of one scan area boundry in degrees; 

Longitude of one scan area boundry In degrees; 

Latitude of other scan area boundry In degrees; 

Longitude of other scan area boundry In degrees; 

Earth central angle of annular radius in degrees; 

Logical variable, TRUE for Inner annulus, FALSE for outer 
annul as search; 

Depth of search area In degrees 


OUTPUT ARGUMENTS: TOTTRF 


TOTTRF Total number of background aircraft traffic found in the 


scan area. 





METHOD: 


A series of ten degree wide subregions Is defined within the scan 
area. The midpoints of these subregions and their areas are used to 
define the density blocks for routine B6T. The final subregion may be 
less than ten degrees In width, being the remainder of the region. The 
area of a subregion Is defined to be depth times the width times the 
degree-to-nautical-mile conversion squared. 

REFERENCES: 

BGTSCN Is called by BRRIER Subroutine. 

BGTSCN calls BGT, BRNG, DGLIM, FLALOB and TRI. 






3.3.4 SUBROUTINE BRRIER 


PURPOSE: 


Determines If a satellite detection annulus edge intersects an 
aircraft barrier. If there Is an intersection, BRRIER tests to see if any 
aircraft are detected In the search swath along the annulus edge and between 
the two points of Intersection. Because aircraft barriers can be annulur, a 
satellite detection annulus edge can Intersect a barrier twice. 


OUTER EDGE 
OF BARRIER 


AIRCRAFT BARRIER 

EDGE OF 
DETECTION 
ANNULUS 


SEARCH SWA 



BARRIER CENTER 


INNER EDGE OF BARRIER 


Aircraft barriers are searched at the inner and outer edge of the annulus. 


INPUTS: 


TARGTS 


BELEVO 


BELEVI 


TPOS (1,NTAR) 
TPOS (2,NTAR) 
TPOS (4,NTAR) 


Array: See BRBLD for description; 

Array: See BRDLD for description; 

Aircraft barrier ID; 

Outer annulus edge at wheel barriers are searched 
(degree earth grazing angle); 

Inner edge of annulus at which aircraft barriers 
are searched (degrees, earth grazing angle); 
Aircraft latitude; 

Aircraft longitude; 

Time since aircraft began course; 

Beamwidth of radar beam; 


tl 







OUTPUTS: 


NACS 

Number of aircraft barriers being searched by 

NTSCH 

current spacecraft; 

Number of aircraft and ship barriers 

being 

Locations of annulus barrier 

searched by current spacecraft; « 

Intersections defining search swaths 

are 

stored In: 



PT (1, NTSCH, NSC) 

Latitude of first Intersection; 


PT (2, NTSCH, NSC) 

Longitude of first Intersection; 


PT (3, NTSCH, NSC) 

Latitude of second Intersection; 


PT (4, NTSCH, NSC) 

Longitude of second Intersection; 


TPOS (1, BROW) 

Latitude of first Intersection; 


TPOS (2, BROW) 

Longitude of first Intersection; 


TPOS (3, BROW) 

Latitude of second Intersection; 


TPOS (5, BROW) 

Longitude of second Intersection; 



Where BROW Is 

- first row In TARGTS array with the barrier ID (for the first swath) 

- second row In TARGTS array with this barrier ID (for the second 
swath) 


/?\ 


Note: TPOS Is changed every time the sensor loop Is Incremented while PT Is 
changed only when the time loop Is Incremented. PT Is used to write 
latitudes and longitudes of swath limits on the SCPLT file. 


RES (1, NISCH) 

RES (2, NISCH) 
RES (3, NISCH) 
RES (4, NISCH) 


Width of search swath, measured In spacecraft 
centered degrees; 

Search cycle time (minutes); 

Available dwell time for a single bean (minutes); 
Width of search swath, measured In rotational 
degrees about suborbital point; 
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run. * vt *&***%***»,»' I 


RES (5, NISCH) 


RES (6, NISCH) 


Annulus edge bordering search swath (earth 
central angle In degrees); 

Number of aircraft (Includes background traffic 
plus threat aircraft) that are in search swath 
and that need to be tracked. 














METHOD: 


BRRIER performs three functions for each aircraft barrier for each 
spacecraft sensor at each time step: 

1. Determines If the Inner or outer edge of the annulus Intersects the 
barrier; thus creating a search swath; 

2. Completes the entries In the RES array for each search swath.; 

3. Determines If there are any threat aircraft In the search swath 
being detected for the first time; 

If there are two search swaths, then the ID for the first one Is the 
odd number Input as the barrier ID In the TARGTS array and the ID of the 
second one Is the Input ID number plus one. 

To determine If there Is a search swath, BRRIER uses CONES to test for 
the Intersections of the annulus edge with the following four circles: 

1. The circle defining the Inner edge of the barrier; 

2. The clrle defining the outer edge of the barrier; 

3. The great circle defining each end of the barrier. See description 
of BRBLD. 

If these Intersections exist, they are then tested against the barrier 
specifications input on the TARGTS array. This Is done by transforming the 
Intersections to a new coordinate system whose pole Is the barrier center and 
whose zero meridian Is the barrier's bisector (or axis of symmetry). After 
this transformation, the colatitude and longitude of each Intersection are the 
distance from the barrier center and the bearing relative to the barrier 
bisector. Then, for an Intersection to be on the barrier, the colatitude must 
be on or between the Inner and outer edges of the barrier and the bearing must 
be less than or equal to the half-angle of the barrier width (see BRBLD). 
Intersections passing these tests are Input to the SRT array for scheduling by 
calling FSTEP. 

The entries In RES are then computed for each swath. Search swath 
width In rotational degrees Is computed by converting the points of 
intersection to a new coordinate system whose pole Is the spacecraft 
suborbital point and whose zero meridian Is along the spacecraft's Inertial 
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velocity vector. The depth of the search swath is the depth of the radar beam 
footprint or the distance traveled by the suborbital point in one simulation 
time step, whichever is larger. The period or time available to scan (or 
search) the swath is the time it takes the annulus edge to travel over a 
distance equal to the depth of the swath. This period is also called the 
search cycle time; to maintain a 10W probability of sweeping over an aircraft 
as it goes through the search swath, the search swath must be scanned at least 
once per cycle time. Available beam dwell time is the ratio of cycle time to 
swath width (In spacecraft centered degrees) times the beamwidth. The number 
of aircraft tracked Is the number of threat aircraft detected (described 
below) plus one-half the number of background aircraft in the swath (because 
It Is assumed that half the aircraft in a barrier are headed towards the 

Soviet side of the barrier and half are headed away from the Soviet side). 

Finally, the tests to determine if any threat aircraft, which are 
Input In the TARGTS array, are detected in the swath are: 

1. Has the aircraft begun its course? If yes, then; 

2. Has the aircraft ended its course? If no, then; 

3. Is the aircraft inside the search swath? If yes, then; 

4. Does the aircraft have a doppler return greater than the minimum 
discernable velocity? 

If all these tests are passed, then the aircraft's position relative 
to the spacecraft is input to the SRT array by calling FSTEP. 






' 3.3.5 SUBROUTINE BRRIRS 

* 

PURPOSE: 

Models detections of ship barriers by radar satellites. A ship search 
swath Is defined by the segment of the outer edge of an annulus that Is Inside 
a ship barrier. Because ship barriers are always circular, there Is never 
more than one search swath per barrier. 

INPUTS: 

SHELEV Annulus edge at which ship searches are conducted. 

Expressed In terms of earth grazing angle (degrees); 
IBR Number of ship barrier being tested; 

BREFS ARRAY See BRBLD for description; 

NSC Number of current spacecraft; 

OUTPUTS: 

NSHS Number of ship barriers being searched by current 

spacecraft; 

NTSCH Number of aircraft and ship barriers being searched by 

current spacecraft; 

Each ship search swath Is bounded by the Intersections of the search annulus 
and the circle describing the barrier. The latitude and longitude (In 
degrees) of each Intersection Is stored as follows: 

TPOS (1, NTAR) Latitude of first Intersection; 

TPOS (2, NTAR) Longitude of first Intersection; 

TPOS (3, NTAR) Latitude of second Intersection; 

TPOS (5, NTAR) Longitude of second Intersection; 

PT (1, NTSCH, NSC) Latitude of first Intersection; 

PT (2, NTSCH, NSC) Longitude of first Intersection; 

PT (3, NTSCH, NSC) Latitude of second intersection; 

PT (4, NTSCH, NSC) Longitude of second Intersection; 






RES (1, NTSCH) 


Width of search swath, measured In spacecraft 
centered degrees; 

RES (2, NTSCH) Search cycle time; 

RES (3, NTSCH) Available dwell time for a single beam; 

RES (4, NTSCH) Width of search swath, measured In rotational 

degrees about the suborbital point; 

RES (5, NTSCH) Annulus edge bordering search swath; expressed as 

an earth central angle (In degrees) from the 
suborbital point; 

RES (6, NTSCH) Not used; 

ISCAN (NTAR, NSC) Set to one If barrier is searched by spacecraft; 
is used on the SCAN output file; 

STBER 1, STBER 1 Bearings (in degrees) from suborbital point to 
Intersections of annulus edge and ship barrier; 

ST Earth central angle between suborbital point and 

annulus edge at which barrier search Is 
performed; 

METHOD: 

Subroutine CIRSCT Is called to compute location of two Intersections 
between annulus edge at which ship barriers are searched and circular ship 
barrier. If there Is no Intersection then CIRSCT sets IPOINT to zero, which 
causes BARRIRS to return Immediately. If the two circles do Intersect, then 
o NSHS and NTSCHS are incremented; 
o The latitude and longitude of the Intersections are 
stored In TPOS and PT; 

o ISCAN (NTAR, NSC) Is set to one, when NTAR Is the 
row of the TAR6TS array describing this barrier 
and NSC Is the number of the current spacecraft; 
o Subroutine FSTEP Is called to enter the two points 
of Intersection into the SRT array for clustering 
and scheduling; 

o The first five entries In RES are computed; 







3.3.6 SUBROUTINE CLUSTR 


PURPOSE: 

Partition all the search areas in the annulus of a single active 
satellite Into groups (or clusters) which can be mechanically scanned by 
moving the antenna feed of the satellite boresight. There are two types 
of search areas, barrier search areas and aircraft reacquisition and 
tracking areas. A barrier search area Is a swath, one beam footprint 
deep, along the edge of the annulus which Is Inside a barrier. Aircraft 
reacquisition and tracking areas are much smaller areas of uncertainty 
which are searched to detect aircraft that have been previously 
detected. For each cluster, the following are computed: 

o boresight cumpolnt 
o background traffic In search areas 
o energy used to scan cluster 


INPUTS: 

CLSWID degrees Angle, centred at suborbital point and 

expressed In rotation degrees about that 
point, defining the maximum width of a 
cluster, l.e., the maximum width that 
can be scanned from a fixed boresight 
position by moving the antenna feed. 

NTD Total number of targets detectble In 

annulus. In the context, targets 
Include: 

o Intersections of barrier 
boundaries with annulus 
o aircraft Inside annulus 
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AA (35,35) A square symmetric zero diagonal array 

containing the suborbital point centered 
angle between all pairs of tarqets 
Inside the annulus. For barrier search 
areas, aircraft not yet detected, and 
ships in barrier search areas, the 
diagonal element is reset from zero to 
the ID number of the search area. 

NSC Number of current spacecraft 

SRT (3,35) Stores for each target in annulus: 

1. Row of target on TARGTS array 

2. Bearing, relative to spacecraft 
velocity bearing, from suborbital 
point to target (degrees) 

3. Spacecraft look anqle, from nadir 
(degrees). 


OUTPUTS: 


CLM(36,10) These arrays can store up to 41 items 

CLM2(5,10) of information about each cluster: 

CLM3(25,10) 

CLM(l.ICL) Cluster number 

CLK(2,ICL) Number of spacecraft viewing cluster 

CLM(3,ICL) Cluster number for this spacecraft 

CLM(48t5,ITCL) Bearings from spacecraft suborbital 

point to targets that are gratest angle 
apart; 








CLM{6&7,ITCL) 


CLM(8,ITCL) 

CLM(9,ITCL) 

CLM(IO.ITCL) 

CLH(11,ITCL) 
CLM(12-36,ITCL) 


Spacecraft look angles to see targets 
that are greatest angle apart. 4-6 are 
for the target tha is closest to the 
suborbital point; 5-7 are for the tarqet 
that is furthest from the suborbital 
point; 

Bearing to boresight AIM point 
Spacecraft look angle to boresight AIM 
point; 

Maximum angular separation of targets in 
cluster; 

Number of targets In cluster 
Number of each target in cluster; 
Initially the number is the target's 
rank in the SRT array. Maximum number 
of targets currently allowed in a 
cluster is 25.; 


CLM2(1,JCL) 

hours 

Earliest time of last detection across 

all targets in the cluster; 

CLM2(2,JCL) 

hours 

Minimum duration of current detection 

Interval across all targets in cluster; 

CLM3(1,ITCL) 


Barrier number; 

CLM3(2,ITCL) 


Earth central angle of search (inner or 
outer edge of annulus) in degrees; 

CLM3(3,ITCL) 


Search cycle time (minutes) from 
barrier; 

CLM3(4,ITCL) 


Search amplitude in spacecraft degrees; 

CLM3(5,ITCL) 


Search time (minutes), hopefully less 
than cycle time; 

CLM3(6,ITCL) 


Search energy (KJ) expended during time 
step; 

CLM3(7,ITCL) 


Background traffic In search area (Inc 
targets); 

CLM3(8,ITCL) 


Time P/l on during search (hrs); 


£ 


€ 






ill lot"' 


-X. 
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CLM3(9,ITCL) 


Number of targets to be reacquired; 

CLM3(10,ITCL) 


Energy (KJ) used to reacquire targets; 

CLM3(11,ITCL) 


Time (minues) to reacquire targets; 

CLM3(12,ITCL) 


Background traffic in reacquisition 
areas (inc targets); 

CLM3(13,ITCL) 


Background traffic in tracking areas 
(inc targets); 

CLM3(14,ITCL) 


Number of raid targets to be tracked; 

CLM3(15,ITCL) 


Track time/cycle time ratio; 

CLM3(16,ITCL) 


Maximum number of targets that can be 

tracked in the available track time; 

CLM3(17,ITCL) 


Track energy (KJ); 

CLM3(18,ITCL) 


Total background traffic encountered; 

CLM3(19,ITCL) 


Time P/L on during tracking (hrs); 

CLM3(20,ITCL) 


Time P/L on during reacquisition (hrs); 

CLM3(21,ITCL) 


Total energy used (KJ) to search, 
reacquire, and track targets; 

CLM3(22,ITCL) 


Total time P/L on during period; 

CLM3(23,ITCL) 


Not used; 

CLM3(24,ITCL) 


Not used; 

CLM3(25,ITCL) 


Not used; 

CLCTR(2,5,8) 

CLCTR(1,IC2,NSC) 

degrees 

Longitude of boresight aimoint; 

CLCTR(2,ICL,NSC) 

degrees 

Latitude of boresight aimpoint; 

CLBG(5,5,9) 

degrees 

Stores spacecraft centered angles 
between all pairs of boresight aimpoints 
in annulus; 

The diagonal, which otherwise would be 
zero, stores: 

CLBG(1,1,NSC) 


Is # if there is a scheduling conflict, 
i.e., if this spacecraft can see more 
than one cluster; 

CLBG(2,2,NSC) 


Number of clusters in annulus; 
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CLBG(3,3,NSC) Number of clusters for which this 

spacecraft is responsible. (Another 
spacecraft may be able, for example, to 
see one of the clusters visible to this 
spacecraft.) 

Total angle to visit all clusters. 


CLBG(4,4,NSC) 







METHOD: 


Subroutine CLUSTR partitions th targets In the AA array Into 
clusters, which are defined and stored In the CLM, CLM2, and CLM3 
arrays. Once all the targets In a cluster are Identified, the boreslght 
almpolnt (In the center of the cluster) and the energy required to scan 
the cluster are computed. 

The following rules guide the building of a cluster: 


1. A cluster contains no more than one barrier search area swath. 

2. The entire swath (defined as the annulus edge between two 
annulus - barrier Intersections) must be In one cluster, even 
If the resulting cluster width is greater than CLSWID. (If 
this occurs, the spacecraft must slew while the antenna feed 
Is moving). 

3. A cluster will generally contain an aircraft barrier or a ship 
barrier and several aircraft track or reacquisition search 
areas which are approximated by the current, true position of 
the aircraft. 

4. Clusters are first built around barrier search swaths and then 
around unclustered aircraft. The first member of a cluster Is 
called the pivot. 

5. If the pivot Is a barrier search swath intersection or as yet 
undetected aircraft or ship appearing in a barrier search 
swath, then all other members of that barrier search swath are 
included In the cluster. CLSWID Is replaced by the width of 
the barrier search swath If it Is greater than the search 
swath. 

6. After the pivot target(s) are in the cluster, all other 
aircraft in the AA array are tested to see If they can also be 
in the cluster. To be in the cluster, the candidate aircraft 
must be less than CLSWID degrees away from all targets already 
In the cluster. 

7. The process of Identifying a pivot and testing all other 
unassigned AA array entries against that pivot continues until 
all entries of the AA array are clustered. 

8. The diagonal element In the AA array for an Individual target 
Is set to 99 when that target has been Included In a cluster. 
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The boreslght Is aimed at a spacecraft look angle that bisects the 
spacecraft centered angle from the Inner to the outer edge of the annulus 
and at bearing angle that bisects the width of the cluster. 

The logic to compute background traffic In search areas and the 
time and energy to scan a cluster Is still developmental and Is not 
correct as of this writing. 
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3.3.7 SUBROUTINE CLSTPP 


PURPOSE: r 

Subroutine CLSTPP (cluster port processor) Is the scheduler for active 
satellites. If an active satellite can view more than one cluster (see 
description of subroutine CLUSTR), then the schedule decides at which cluster 
to aim the boreslght. 


INPUTS: 

CLBG (5,5,9) 


CLM (36,10) 


See description of subroutine CLUSTR. CLSTPP 
checks to see If two satellites see the same 


cluster. In which case the third diagonal element 
Is decremented for one of the spacecraft and the 
first diagonal element may be set to a blantk. 

See description in CLUSTR. For each cluster 
(l.e., for each row on the CLM array), CLSTPP 

uses: 


CLM (l.ICL): 

Cluster number; 

CLM (2,ICL): 

Number of spacecraft 
viewing cluster; 

CLM (3,ICL): 

Cluster number for this 

spacecraft; 

CLM (11,ICL): 

Number of targets in 
this cluster; 

CLM (12-36): 

Target number (l.e., row 
on targets array) of 
each target In cluster; 


Note: In this context, "target" can mean the endpoints of a barrier search 
swath. 


CLM3 (1,ICL) 


Earliest time of last detection across all 
targets In cluster. 





CLM2 (2,ICL) 


CLM3 (1,ICL) 
XLN (NTAR.l) 
XLN (NTAR.U) 


NSC I 

OUTPUTS: 

SCSCH (9,8) 

SCSCH (l.NSL) 

SCSCH (2,NSC) 

SCSCH (3,NSC) 

SCSCH (4,NSC) 

SCSCH (5,NSC) 

SCSCH (6 and 7.NSC) 

SCSCH (8,NSC) 

SCSCH (9,NSC) 

ICLST (9, 10) 


Minimum duration of current detection Interval 
across all targets In cluster. 

10 number of barrier search swath 
Time (hours) at which last detection occurred. 
Time (hours) at which current detection began; Is 
set to zero If target was not detected at last 
time step. 

Total number of clusters seen by all active 
sensors (e.g. active spacecraft, 0T4 B radars, 
and DEW line) at current time step. (By 
definition, there Is no more than one cluster per 
fixed sensor.) 

The number of active sensors. 


Shows how each active sensor Is scheduled. 

1 If sensor sees no more than one cluster; 

0 If sensor sees none or several clusters; 

1 If sensor has been assigned to a cluster; 

0 If sensor has not been assigned to a cluster; 
Row In CLM/CLM2/CLM3 arrays containing cluster to 
which this sensor has been assigned; 

Number of cluster to which spacecraft has been 
assigned. This number can be less than SCSCH 
(3,NSC), but only If two or more sensors see the 
same cluster (l.e., have matching clusters). 
Number of targets In cluster; 

First and last rows In CLM/CLM2/CLM3 containing 
clusters that can be seen by this spacecraft; 
Earliest, time of last detection across all 
targets In assigned cluster; 

Minimum duration of current detection Interval 
across all targets In assigned cluster; 

Shows which active sensors are assigned to each 
cluster; 
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ICLST (1,ICL) Total number of active sensors assigned to this 
cluster; 

ICLST (2-9,ICL) Active sensor ID numbers of sensors scheduled to 
view this cluster; 

ICLWID (8) Width of cluster (In terms of rotation degrees 

about spacecraft suborbital point) which has been 
scheduled for each spacecraft; 

METHOO: 

CLSTPP assigns or schedules active sensors to look at one cluster per 
time step. For OTH B radar and DEW Line this Is trivial, because by 
definition OTH B and the DEW Line each see only one cluster. Scheduling 
spacecraft Is harder, because. In general, each spacecraft has at each time 
step more than one cluster (In Its annulus) that needs to be scanned. Recall 
that a cluster Is a group of aircraft and up to one aircraft or ship barrier 
search swath that can be scanned from a fixed boreslght position by moving the 
antenna feed. When an active spacecraft has a scheduling conflict, the 
conflict Is resolved based on the following priorities (listed In descending 
order of Importance). 

1. Do not sche&le a cluster If It Is already scheduled to another 
sensor (Implies that two sensor see the same cluster) unless the 
current sensor has In Its cluster a target which: 

- Is not seen by the other sensor; 

- has longer detection gap than the targets by the other sensor; 

- has a shorter current detection period than that of any of the 
targets seen by the other cluster; 

2. Choose the cluster that has target with the longests detection gap 
unless the maximum gap Is less than .9 hrs. 

3. Choose the cluster which has the target with the shortest detection 
period If: 

- the target was observed at the last time step; 

> The detection period Is less than .3 hrs; 




4. Choose the cluster that has the most targets. 

It occasionally happens that a sensor is not scheduled to see a 
cluster - all the clusters it sees are already scheduled to other sensors. 

The scheduler assigns sensors In ascending order -1, 2, ...» NSCI; it does not 
unassign nor reassign sensors. (This would be a desirable Improvement.) When 
a sensor Is scheduled to a particular cluster, this Is recorded In the SCSCH 
array. 


3.3.8 Subroutine ECLPSF 
PURPOSE: 

This routine Is responsible for updating the SUNVEC data to the 
current simulation date/tlme. It Is called within the early stages of the 
Time Loop of COV. The SUNVEC Is used to determine the eclipse status of 
the spacecraft within that Time Loop. 

INPUT ARGUMENTS: 

DATE The current Julian date and tine, the latter expressed as 
a fraction of a day. 

OUTPUT ARGUMENTS: 

SUNVEC The x-, y- and z- unit vector for the sun's position 
relative to the Earth. 


METHOD: 

The algorithm Is embeded In SUN and the other utility routines 

called. 

REFERENCES: 

ECLPSF Is called by COV. It calls the following utility 
routines: SUN, OTOC, EQEC, and RMMUL. 
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3.3.9 Function ECLPSF 
PURPOSE: 

ECLPSF Is a logical function Intended to return a TRUE value Mhen 
the spacecraft position Input Is In solar eclipse during the current time, 
and otherwise, FALSE. 

INPUT ARGUMENTS: 

SUNVEC The x-, y- and z- unit vector of the sun's position 
relative to the Earth at the current time. 

SCLON Longitude of suborbital point in degrees 
SCLAT Latitude of suborbital point In degrees 
SCRAO Orbital radius In feet. 

OUTPUT ARGUMENTS: 

ECLIPSE A logical variable, set TRUE for eclipse, otherwise, 

FALSE. 

METHOD: 

The unit vector for the spacecraft Is computed and the dot product 
of this vector with the solar unit vector Is taken. From the comparison 
of the result with the horizon of the spacecraft, the eclipse state can be 
determined. 

REFERENCES: 

The routine Is used in COV. It calls DOT and MAGN. 





3.3.10 SUBROUTINE FSTEP 


PURPOSE: 

Performs the first step of active satellite clustering and scheduling, 
which is to store the bearing and elevation angles of detectable search areas 
in the SRT array in order of their bearing relative to the spacecraft's 
velocity heading. 

INPUTS: (all angles In degrees) 

Earth central angle from suborbital point to a search 
area location; 

Bearing (from north) from suborbital point to search 
area location; 

Spacecraft velocity bearing; 

Search area ID number (row on TARGTS array that 
describes the search area). A search area can be a 
swath through an aircraft or ship barrier or the 
estimated current position of an aircraft. 

OUTPUTS: (all angles In degrees) 

SRT(1, I) Search area ID number; 

SRT(2, I) Bearing, relative to spacecraft velocity heading, from 

spacecraft suborbital point to search area location; 
SRT(3, I) Spacecraft look angle from nadir to target; 

Note: all detectable search areas are entered In the SRT array In ascending 
order based on bearing from suborbital point to search area. Bearings can 
vary from -180° to + 180°. 

NTO Number of search areas detectable by current 

spacecraft; 


ST 

STBER 

SVBER 

NTAR 




! 


METHOO: 

This subroutine Is an intermediate step which is prerequisite for 
subroutine CLUSTR. Subroutine CLUSTR groups detectable search areas for an 
active satellite into groups (called clusters) which can be scanned from a 
fixed boresight position by moving the antenna feed. To do this, CLUSTR 
requires that all detectable search areas be stored In clockwise order around 
the annulus. FSTEP puts search areas in this order. A search area is: 

o An aircraft barrier search swath; 

o A ship barrier search swath; 

o The current position of an aircraft that has already 
been initially detected; 

For each search swath, two points are stored in SRT, namely, the two 
intersections of the barrier with the annulus edge bordering the search 
swath. Currently, the model approximates a search area to reacquire a 
previously detected aircraft by the aircraft's current, actual position. This 
approximation is not valid if the aircraft maneuvers between detections; the 
search area should be based on an a position extrapolated from the last 
observed position and should have some area of uncertainty. 

A search area Is detectable If it Is In the spacecraft's detection 

annulus. 


i 


i 


i 








3.3.11 SUBROUTINE FXSENS 


PURPOSE: £ 

Models detections of aircraft by fixed sensors, e.g., Over-the-Horlzon 
Backscatter (OTH B) radar and the DEW Line. The approach is to model the 
areas covered by these two systems as annular sectors and to compute when 
aircraft are in these sectors. 


INPUTS: 

(All angles in degrees;) 


SLAT 

SLON 

SHO 

SWD 


SOR 


Latitude of annulus center; 

Longitude of annulus center; 

Bearing from north of axis of symmetry (or 
bisector) of annular sector; 

Angulur width of annular sector (half angle); 
Distance from annulus center to Inner edge of 
annular sector (n.ml.); 

01 stance from annulus center to outer edge of 
annulus sector; 


TARGTS (NTAR, 1) Greater than zero if target is an aircraft; 

TARGTS (NTAR, 7) Time at which aircraft begins flight; 


TP0S(1, NTAR) Aircraft's current latitude; 

TP0S(2, NTAR) Aircraft's current longitude; 

TPOS(3, NTAR) Aircraft’s current bearing; 

TP0S(4, NTAR) Time since aircraft started course; 


OUTPUTS: 

CLM(1, J) Cluster number (usually equal to J); 

CLM(2, J) Sensor number; 

CLM(3, J) Number of clusters seen by this sensor (Is always 

one for OTH B radar and DEW Line); 

Number of aircraft seen by this sensor; 


CLM(11, J) 




Aircraft ID (l.e., row in TARGTS array) for each 
aircraft seen by sensor; 


To determine If an aircraft is in the coverage area of an OTH B radara 
or the DEW Line, FXSENS checks all rows in the TARGTS array to find if: 

o row represents an aircraft; if yes, then; 
o has the aircraft begun its flight? if yes, then; 
o has the aircraft ended its flight; if not, then the 
current position of the aircraft Is transformed 
to a new coordinate system whose pole is the 
center of the annular coverage sector (e.g. the 
OTH B radar transmitter Is the center of the 
annular sector covered by the OTH B radar system) 
and whose zero meridian Is the axis of symmetry 
of the annular sector. In the new coordinate 
system, the aircraft's colatitude is its distance 
■% from the center and its longitude is Its angulur 

separation from the center of the annular 
coverage sector. The tests to determine if the 
aircraft is in the annular coverage sector are, 
therefore: 

o Is aircraft's new colatitude greater than the 
distance to the inner edge of the coverage 
sector; If yes, then; 

o Is the new colatitude less than the distance to the 
outer edge of the coverage sector; If yes, then; 
o Is the absolute value of the aircraft's new 

longitude less than the haif-angle width of the 
coverage sector; 

If all these tests are passed, then the aircraft Is In the coverage 
area of the current fixed sensor. 

If the fixed sensor sees one or more aircraft, then FXSENS creates a 
cluster, similar to the clusters created for spacecraft by subroutine 
CLUSTR.. However, whereas for spacecraft there may be several clusters of 


CLM (12-35, J) 


ME^Hof: 
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which only one can be observed at a time; for OTH B radar or the DEW Line, 
there is only one cluster, i.e., all aircraft in the coverage area can be 
seen in one simulator time step. A cluster is created by making the entries 
shown under outputs in a new row of the CLM array. This makes it unnecessary 
to call CLUSTR for fixed sensors. However, subroutine CLSTPP, which schedules 
spacecraft to look at the highest priority cluster, uses the CLM array to 
assign spacecraft to clusters which do not have aircraft that are seen by a 
fixed sensor. Thus, fixed sensors can simplify spacecraft scheduling. 


# 
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3.3.12 SUBROUTINE JMEQ 


PURPOSE: 

Computes an aircraft's velocity component along the line of sight from 
the aircraft tothe spacecraft. 

INPUTS: (all angles in degrees) 

ST Earth central angle from suborbital point to aircraft 

position; 

DBER Suborbital-point-centered angle between spacecraft 

velocity bearing and great circle to aircraft 
position; 

VT Aircraft velocity (knots) 

OUTPUTS: 

DDOP Aircraft's velocity component along the line of sight 

from the aircraft to the spacecraft (knots); 

METHOD: 

Compute the aircraft centered angle between the aircraft velocity 
vector and the aircraft-to-spacecraft line-of-sight radial. This is done 
using a right spherical triangle on a unit sphere whose center Is the 
aircraft. Two adjacent (and perpendicular) sides are the arcs corresponding 
to DBER and the elevation angle (above the aircraft's local horizon) to the 
spacecraft. The third arc length corresponds to the angle between the 
velocity vector and the line-of-siqht from the aircraft to the spacecraft. 

The aircraft velocity component along the line-of-sight is Its speed times the 
cosine of the angle just computed. 
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3.3.13 


SUBROUTINE PWRBAL 


PURPOSE: 

This routine Is responsible for maintaining the state of charge of 
the battery. The state of charge is updated each time increment based 
upon the amount of time the payload Is used during the period and the 
solar eclipse state. 

INPUT ARGUMENTS: 

The following parameters are input via the calling sequence: 

TUSED Amount of time payload was used in hours. 

ECLIPS Logical variable indicating eclipse state, TRUE for 
eclipse; FALSE, otherwise. 

SOC State of charge, as a fraction, at start of time 
Interval. 


The following parameters are input via Common X2 or X8 


BATCAP Maximum battery capacity in kilowatt hours. 
r SOLAR Power input from solar panels in kilowatts 
PPLON Power consumed while payload is on, including power 
consumption of platform, in kilowatts. 

PPLOFF Power consumed while payload Is off. Including power 
consumption of platform, in kilowatts. 

DUTCYC Duty cycle of payload. 

OTIME Time interval In hours. 

OUTPUT ARGUMENTS: 

SOC The state of charge at the end of the time interval 







METHOD: 

The input power Is determined by the eclipse state, 0 for TRUE, 
PSOLAR for FALSE. If the payload was used during the period, the 
consumption rate is set to PPLON for a period of DTYCYC * TUSED. The 
period the payload was off is either the entire Interval or the remainder 
of the interval (DTIME-DTYCYC * TUSED) and power consumption, PPLOFF. The 
net power is then computed (Power in minus Power out) as well as an ETA 
factor. If the battery is discharging, the change in state of charge is 
the net change times the period divided by ETA times battery capacity. If 
the battery Is charging, the change In state of charge is equal to net 
change times ETA times .7576 times the period divided by the battery 
capacity. The new SOC Is then the old SOC plus the change in state of 
charge. If the payload was used, there Is a period that it was not used 
during the interval. The power consumption during that period is computed 
next as already described for the payload off. 

REFERENCES: 

This routine is called by COV. 
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3.4 OUTPUT 


3.4.1 SUBROUTINE BGNCRS 
PURPOSE: 

BGNCRS calls SETL EPSl before moving target begins course so that 
coverage statistics computed by subroutine summary and the post-processor 
STATCV Include gaps and detections occurring at the beginning of the 
target's course. (Note that. In the statistics, the Initial gap or 
detection will be EPSl plus the duration of the gap or detection.) 

ENDCRS calls SETL EPSl after moving target begins course so that 
coverage statistics computed by subroutine summary and the post-processor 
STATCV Include gaps and detections occurring at the end of the target's 
course. (Note that. In the statistics,the final gap or detection will be 
EPSl dIus the duration of the gap or detection.) 

INPUTS: 


STPROF 

hours 

Time (since start of simulation) that moving 
target begins course. 

EDPROF 

hours 

Duration of target's course. 

EPSl 

hours 

See above. Normally stt to .0001 hours. 

NSMCON 

- 

Number of target-ground station sets on 
ISMTAB array. 

NTRGTT 

- 

Number of targets on the TARGTS/ITRGTS array 
plus NSMCON. 

DTIME 

hours 

Time steo size for simulation. 

TIME 

hours 

Current time being simulated. 

OUTPUTS: 




Added to TIME; Is the time at which SETL Is 
called and results In SETL being called EPSl 
before course begins or EPSl after course 
ends. Is reset to 999 at end of subroutine. 


TAOJ 


hours 











METHOD: 

The main program, COV, calls BGNCRS/ENDCRS two time steps 
before/after course begins/ends. For 8GNCRS, a modular divide determines 
how much time elapses in the next time step before the course begins; EPS1 
Is subtracted from the time period to compute TADJ. For ENDCRS, a modular 
divide determines how much time elapsed In the previous time step before 
the course ended; EPS1 Is added to this period to compute TADJ. After 
SETL Is called, TADJ Is reset to 999. 

REFERENCES: 


Calling routines: COV 
Called routines: SETL 






3.4.2 Subroutine CTOALF 


PURPOSE: 

This routine is responsible for converting the numeric cluster 
number found in ISCAN array to either an alpha or a digit (in alpha 
form). The SCSCH(4, I) volumn Is also converted to alpha. The ISA vector 
Is set from the rounded value of the CLBG(5,5,I) vector. These values are 
used in the SCAN printout. 

INPUT ARGUMENTS: 

The following Inputs are passed via the calling sequence: 

NSC1 Number of aopropriate S/C; 

NTRGTS Number of targets; 

SCSCH Array whose (4, I) column represents the selected 
cluster id; 

The ISCAN array and the CLBG array are passed via the X4 Common. 
OUTPUT ARGUMENTS: 

The following outputs are passed via the calling sequence: 

ISCSCH Array whose (4, I) column represents the converted ID for 

the selected cluster; 

ISA Vector containing the rounded value of the CLBG (5, 5, I) 
column. 

The ISCAN array is also used for output. 


METHOD: 


ISA Is derived from the following: 

ISA(I) - IFIX(.5 + CLBG(5, 5, I)) 

ISCAH Is derived from the following: 

If ISCAN (J,I) - 0 

-ISCAN(J,I) as an alpha coded digit (1-9) 
ISCAN(J.I) If ISCAN(J,I) < 0 

•A'-'Z' If 0 < ISCAN(J.I) < 26 
If ISCAN(J,I) > 26 

ISCSCH(4,I) Is transformed In the same way ISCAN Is. 
REFERENCES: 

This routine Is called from COV. 


.1 
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3.4.3 SUBROUTINE SETL 


PURPOSE: 

Gathers statistics on the duration of detections and of gaps between 
detections. At the start of each new detection, SETL also writes on unit 8 
the time at which the detection begins and the duration of the coverage gap 
since the last detection. 

INPUTS: 

\ 

TfNE Current clock time (hours); 

DTIM^ Simulator time step (hours); 

DTIMX Larger time step if two types of spacecraft are 

modeled; 

SATTAK Time at which current replication started (hours); 

TADJ Amount by which time has been adjusted by ENCRS or 

BGNCRS; 

NTAR I (row on TARGTS array) for object being detected; 

IPL0T2 Must be equal to or greater than one to cause a write 

on unit 8; 

OELTIM Time since aircraft began flight or since start of 

replication (for fixed targets); 

NUMCT Number of simultaneous contacts required for 

detection; 

OUTPUTS: (all times In hours) 

XLN(I, 1) Time at which last contact occurred; 

XLN(I, 2) Maximum time out of contact; 

XLN(I, 3) End time of maximum time out of contact; 

XLN(I, 4) Maximum time out of contact; 

XLN(I, 5) End of minimum time out of contact; 

XLN(I, 6) Total time out of contact; 

XLN(I, 7) Number of contacts; 

XLN(I, 8) Sum of squares of down times; 

XLN(I, 9) Time at which first detection began; 










XLN(I, 10) Time at which last detection began; 

XLN(I, 11) Tine at which current contact began (zero If currently 

out of contact); 

Data written on unit 8: 

NTAR 
IATTAK 
DLTIM 
DT 

ICLS (NTAR) 

METHOD: 

SETL Is called everytlme a ground station, an aircraft, or a barrier 
Is detected. If the detection Is the first detection after a gap, then the 
SLM array Is updated and unit 8 written on; otherwise a RETURN Is executed. 
Special tests are required If a detection occurs exactly at the start of a 
replication or at the first time step Into a replication. BGNCRS and ENDCRS 
call SETL just before and after aircraft flights so gaps before the first 
detection and after the last detection are counted In the statistics. 

However, when ENDCRS calls SETL, DT Is always less than the time step, which 
would cause an Immediate return. However, this return Is prevented If TADJ Is 
negative, and TADJ Is negative only If SETL Is called from ENDCRS. 

Special format statements are used to write a zero before single digit 
positive numbers so sort utilities will sort on a zero Instead of a blank. 


ID of aircraft or barrier being detected; 
Replication number; 

DELTIM (see above); 

Time since last contact ended; 

1 If aircraft Is classified; 0 If otherwise; 







3.4.4 Subroutine STATPT 


PURPOSE: 

This routine calls STHIST for each of five terms to printout the 
post-run statistics. The terms are abbreviated SA, CW, SC, EB, and RC. 

INPUT ARGUMENTS: 

The STATS Common provides all the Inputs to this routine. For each 
of the five terms, the following data are provided where ** Is an 
abbreviation for one of the five terms: 


NCNT** 

Total of all frequency counts 

**MIN 

Minimum value 

**MAX 

Maximum value 

**SUM 

Sum of values 

**SUM2 

Sum of values squared 

I**VEC 

Frequency count vector 


OUTPUT ARGUMENTS: 

For each term, the following data are returned: 

**SUM Mean 

**SUM2 Stand deviation 


METHOD: 

The routine contains data statements for the labels for each of the 
five terms. The routine calls STHIST for each term after heading each 
term's report. 

REFERENCES: 

STATPT Is called by COV. It calls STHIST. 


vN 3 
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3.4.5 Subroutine STATUD 


PURPOSE: 

This routine Is responsible for maintaining the statistics for each 
of the five terms. The count, minimum, maximum, sum and sum square are 
maintained and the frequency bins are Incremented as appropriate. The 
terms are abbreviated SA, CW, SC, EB and RC. 

INPUT ARGUMENTS: 

The calling sequence provides the fallowing arguments: 

MODE The term identifier (1 for SA, 2 for CW, 3 for SC, 4 for 
EB and 5 for RC) 

VALUE The value of the term. 

The STATS Common Is used to hold these statistics, as defined in 

STATUP. 

OUTPUT ARGUMENTS: 

The updated STATS Common Is the only output for this routine. 


METHOD: 

Depending upon the value of MODE, the statistics are updated for 
the term. For the basic statistics, the formulae are as follows where ** 
represents the term's abbreviation: 

NCNT** * NCNT** + 1 
**MIN = AMIN1(**MIN, VALUE) 

**MAX * AMAX1(**MAX, VALUE) 

**SUM « **SUM + VALUE 
**SUM2 » **SUM2 + VALUE * VALUE 
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For this frequency array, the maintenance is dependent upon MODE. 
For MODE = 1, 

I = MIN0(18, MAX0(1, 1 + IFIX(-1.E-10 + VALUE/10.))) 

For MODE = 2, 

I * MIN0(11, MAX0(1, 1 + I FIX(-1.E-10 + VALUE/10.))) 

For MODE = 3, 

I = MINO(12, MAX0(2, 2 + TFIX(-1.E-10 + VALUE/10.))) 
or I = 1 if VALUE < 0 
For MODE * 4, 

I = MIN0(11, MAX0{1, 1 + IFIX(-1.E-10 + VALUE/10.))) 

For MODE = 5, 

I = MIN0(9, MAX0(1, 1 + IFIX(-1.E-10 + VALUE/25.))) 
and I***VEC (I) * I***VEC(I) + 1 

REFERENCES: 


The routine is called by COV. 



3.4.6 SUBROUTINE SMMRY 


PURPOSE: 

Computes statistics summarizing one replication using data stored in 
the XLN array. 

INPUTS 

XLN(I, 1) Time at which last contact occurred; 

XLN(I, 2) Maximum time out of contact; 

XLN(I, 3) End time of maximum time out of contact; 

XLN(I, 4) Maximum time out of contact; 

XLN(I, 5) End of minimum time out of contact; 

XLN(I, 6) Total time out of contact; 

XLN(1, 7) Number of contacts; 

XLN(I, 8) Sum of squares of down times; 

XLN(I, 9) Time at which first detection began; 

XLN(I, 10) Time at which last detection began; 

XLN(I, 11) Time at which current contact began (zero If currently 

out of contact); 

NTR6TS Number of aircraft and barrier Input In TARGTS array; 

OUTPUTS: (all times In hours) 

All outputs are written on unit 6. 

The following Is computed for each target: 

Maximum value of down time; 

Time at which maximum down time ended; 

Minimum value of down time; 

Time at which the minimum down time ended; 

Total down time; 

Number of contacts; 

Mean time between contacts; 

Standard deviation of down times; 
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Percent of time covered; 

Average revisits per day; 

Start time of first contact; 

Start time of last contact; 

In addition the minimum, average, and maximum values across all 
targets are compiled for: 

Maximum time between contacts; 

Minimum time between contacts; 

Average time between contacts; 

Standard devlaton of down time; 

Average duration of contacts; 

Percent coverage; 

Average revisits per day; 

Average revisit interval; 

METHOD: 

Method is implied by the outputs. 
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3.5 UTILITIES 


3.5.1 SUBROUTINE ANGLED 
PURPOSE: 

Computes the angular separation between two points In degrees 

INPUTS: 

A(3), 8(3) — Unit vectors. 

OUTPUT: 

ANG — Angular separation between A and B In degrees. 

METHOD: 

Standard trigonometry. 

REFERENCES: 

No external references 
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3.5.2 SUBROUTINE BRNG 


PURPOSE: Given two pts., A and B, find the gearing from pt. A to pt. B 
In radians and In degrees. 


INPUT ARGUMENTS: 


(.attitude of pt. A (degrees) 
Longitude of pt. A (degrees) 
Lattltude of pt. B (degrees) 
Longitude of pt. B (degrees) 


OUTPUT ARGUMENTS: 


RBRNG Bearing from pt. A to pt. B In radians 

DBRNG Bearing from pt. A to pt. B in degrees 


METHOD: 


Given two pts. A and B, 

1) IF ALAT = 90.0 or BLAT * -90° THEN the bearing Is 
equal to PI. 

2) IF ALAT = -90 of BLAT * 90° THEN the bearing Is 
equal to 0. 

3) Otherwise 

a) Find the distance from pt. A to pt. B. 

b) IF the difference between ALON and BLON Is 
not 180°, then the bearing from pt. A to pt. 
B is defined as follows: 

cos(Bearing) - sin(BLAT) - s1n(ALAT)*cos(DIST ) 
cos(ALAt)*s1n(DI$t) 


DIST Distance from A to B in radians 
Bearing Bearing from pt. A to pt. B 








c) IF the difference between ALON and BLON equals 
180° then the bearing Is defined as follows: 


IF AUT and BLAT are approximately equal then 

If ALAT BUT then bearing - f 
If AUT BUT then bearing ■ 0 

IF. AUT and BUT do not approximate each other then 

If AUT 0. then bearing ■ 0 
If AUT<0 then bearing * 1 

REFERENCES: 

Subroutine BRNG Is used by COVERAGE and PROFILE. 





3.5.3 SUBROUTINE CIRSCT 


PURPOSE: 

Given two points on the earth, each with a different radius, compute 
the points of intersection. 

INPUTS: 

LATA, LONA — Latitude and longitude on point A. 

RADA — Radius of circle at point A in NMI. 

LATB, LONB — Latitude and longitude of point B. 

RADB — Radius of circle at point B in NMI. 

OUTPUTS: 

LATl, L0N1 — Intersect point 1. 

LAT2, L0N2 — Intersect point 2. 

IPOINT — Number of intersections (0 or 2). 


METHOD: 


This routine will never allow only one point of Intersection. If the 
two points are located less than 0.6 degrees of each other, then It Is 
assumed there Is no Intersection. 

REFERENCES: 

The following subroutines are called: 

1. VEC 

2. CONES 

3. ANGLED 

4. RADECM 
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3.5.4 


SUBROUTINE COMBLK 


PURPOSE: 


To set up constants and units for COVERAGE. 


3.5.5 SUBROUTINE CONES 


PURPOSE: 


Computes the two unit vectors defined by the intersection of two given 


INPUTS: 


A(3) — Unit vector, axis of cone 1. 

B(3) — Unit vector, axis of cone 2. 

C0S1 — Cosine of the half angle for cone 1. 
C0S2 — Cosine of the half angle for cone 2. 


OUTPUTS: 

VI(3), V2(3) — The two unit vectors defined by the Intersection of 
the two cones. If the two cones do not Intersect, then Vl(l) Is set to 2.0. 

METHOD: 

Standard trigonometry 
REFERENCES: 


No external references. 
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.5.6 


SUBROUTINE DIS 


PURPOSE: Given two pts., A and B, find the distance from pt. A to pt. B 
In radians, degrees, and nautical miles. 

INPUT ARGUMENTS: 

ALAT Lattltude of pt. A 

ALON Longitude of pt. A 

BLAT Lattltude of pt. B 

IOPT * 1 Implies coordinates of A and B are Input In degrees 

* 2 Implies coordinates of A and B are Input In radians 

OUTPUTS ARGUMENTS: 

RDIST Distance from pt. A to pt. B In radians 

DDIST Distance from pt. A to pt. B In degrees 

DISTNM Distance from pt. A to pt. B In nautical miles 

METHOD: 

If pt. A and pt. B are two points on a unit sphere, then the 
arc-length distance between them Is defined as follows 

cos(DIST) = s1n(ALAT)s1n(BLAT) + cos(ALAT)cos(BLAT)cos(ALON-BLON) 
where 0°<DIST<180° 

output DIST In radians, degrees, and nautical miles. 

REFERENCES: 

Subroutine DIST Is used by COVERAGE and PROFILE. 








3.5.7 SUBROUTINE FLALOB 

PURPOSE: Given the coordinates of a pt. A, the bearing from pt. a to 
another pt., B, find the coordinates of pt. B. 

INPUT ARGUMENTS: 

ALAT Lattitude of pt. A 

ALON Longitude of pt. A 

OIST Distance from pt. A to pt. B 

BRNG Bearing from pt. A to pt. B 

IOPT * 1 Implies that all Inputs and outputs are In degrees 
3 2 Implies that all Inputs and outputs are In radians 

OUTPUT ARGUMENTS: 

BLAT Lattitude of pt. B 

j§ BLON Longitude of pt. B 

METHOD: 

The coordinates of pt. B can be defined as follows: 

sln(BLAT) 3 cos(BRNG)cos(ALAT)s1n(DIST) + s1n(ALAT)cos(DIST) 

BLON * SIGN(DLONG, BRNG) 

where DLONG Is the difference between ALON and BLON 
cos(DLONG) 3 cos (DIST|^s1n(BLATj> sin (ALAT) 

BLON Is then adjusted such that 
0<|BLON|<180° 

REFERENCES: 

FLALOB Is used by COVERAGE. 

© 
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3.5.8 Subroutine HISTGM 


PURPOSE: 

HISTGM is a utility routine responsible for generating a 
horizontally arrayed histogram bar chart on a character printer (as 
opposed to a graphics display device). Both self-scale and user-scale 
options are available. Both percent and absolute levels are printed. 

INPUT ARGUMENTS: 

NCOUNT Total of all frequency counts; 

IFREQ Array of frequency counts (rows); 

NROW Number of bins or rows; 

MAX Scaling term: 

> 0 for the maximum scale (frequency for any row); 

± 0 for self-scale; 

LABEL Array of eight character per row alphanumeric row labels; 
IUNIT The FORTRAN logical unit number of the file on which the 

histogram is to be written. 

OUTPUT ARGUMENTS: 

None 


METHOD: 

If MAX is less than or equal to zero, the routine finds the maximum 
frequency for all the rows. 

All frequencies will be scaled by 50./maximum. 

The histogram is a bar chart whose length is computed by: 

J - TFIX(FLOAT(IFREQ(I))*(50./Maximum) + .5) 

J * MIN0(50, MAX0(0,J)) 
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3.5.9 SUBROUTINE JDORGF 

PURPOSE: 

To convert a date from Julian to year, month, day, hours, minutes, 
seconds or from year, month, day, hours, minutes, seconds to julian date. 


3.5.10 SUBROUTINE LINTRP 
PURPOSE: 

LINTRP is a utility routine responsible for performing a bounded 
linear interpolation of the value of a dependent variable given the 
Independent variable and a table of dependent and Independent variables. 
The Independent variable must be monitopically Increasing. 

INPUT ARGUMENTS: 

XI Independent variable 

X Table of Independent variables 

Y Table of dependent variables 

N Length of tables 

OUTPUT ARGUMENTS: 

Y1 Dependent variable associated with XI. 

METHOD: 

Yo If XI i Xo 

Yl - Y1 +(Y1 - Y1-1 (XI - X1-1)/(X1 - X1-1) for X1-1 < XI < XI 
Yn if XI > Xn 
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3.5.11 


SUBROUTINE MAGN 


PURPOSE: 

To compute the magnitude of vector A. 

INPUTS: A A three element array. 

OUTPUTS: 

B Magnitude of A; 

C Magnitude of A squared; 

METHOD: 

Magnitude of vector A Is the square root of the sum of the squares of 
the three elements of A. 

REFERENCES: 


MAGN Is used by COVERAGE 










3.5.12 SUBROUTINE RADECM 
PURPOSE: 

Computes the right ascension and declination of a given vector, and 
Its magnitude. 

INPUTS: 

ALPHA — Right Ascension In degrees. • 

DELTA — Declination In degrees. 

VMAG — Magnitude of V. 

METHOD: 

Standard trigonometry 
REFERENCES: 

No external references. 
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3.5.13 


SUBROUTINE RMMUL 


■53537*: 










PURPOSE: 

To find the product of two matrices. 

INPUTS: 

A 
B 

OUTPUTS: 

C 


REFERENCES: 

RMMUL is used by COVERAGE. 

C 


Matrix of size L x N; 
Matrix of size M x N; 


Product of matrices A and B; 
(Matrix of size L x N) 
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3.5.14 SUBROUTINE SPHTRN 

PURPOSE: Given the coordinates of a point P In the old system and the 
coordinates of a new pole In the old coordinate system, find 
the coordinates of point P In the new system. 

INPUT ARGUMENTS: 

PHI Longitude of pt. In old system 
TH Colatitude of pt. In old system 

PHIO Longitude of the new pole 
THO Colatitude of the new pole 
PPHIO Bearing of the azimuthal reference 

OUTPUT ARGUMENTS: 

PPHI Longitude of pt. In the new system 
PTH Colatitude of pt. in the new system 

METHOD: 

The new azimuthal reference Is 9 ° 0. Therefore the colatitude 
of the pt. P In the new system can be defined as follows: 

cos(PTH) = cos(TH0)cos(TH0)cos(TH)+s1n(TH0)s1n(TH)cos(PHIO) 
where PTH Is the colatitude of pt. P In the new system. 
sin(PPHI) « -s1n(PHI0)s1n(TH)/s1n(PTH) 

REFERENCES: 

SPHTRN Is used by COVERAGE. 
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3.5.15 Subroutine STHIST 


PURPOSE: 

Tills utility routine Is responsible for computing the mean and 
standard deviation of a set of data and outputting the minimum and 
maximum, total count, mean and standard deviation and frequency histogram. 

INPUT ARGUMENTS: 

NCOUNT Total of all frequency counts, N; 

IFREQ Array of frequency counts; 

NROW Number of rows or bins Into which frequency has been 
divided; 

MAX Explicit scale for histogram If greater than 0; 

LABEL Eight character row labels for histogram; 

IUNIT FORTRAN logical unit number for output; 

XMIN Minimum value of X; 

XMAX Maximum value of X; 

SUM SUM of X; 

SUM2 SUM of X squared; 

OUTPUT ARGUMENTS: 

SUM Mean value of X 

SUM2 Standard deviation of X 

METHOD: 

The mean and standard deviation of X are computed as follows for 
NCOUNT great than 1: 


SUM - SlM/FL0AT(NCOUNT) 

SUM2 - SQRT((SUM2 - FLOAT(NCOUNT)*SUM*SUM)/FLOAT(NCOUNT-1)) 












*7^7^'.v--.'*-V*^ y-r* _. 


The routine writes the NCOUNT, minimum, maximum, mean and standard 
deviation on device IUNIT. It then uses routine HISTGM to generate a 
frequency histogram. 

REFERENCES: 

STHIST Is called by STATPT routine and It calls HISTGM. 
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3.5.16 


SUBROUTINE SUN 


PURPOSE: 

To compute the orbital elements of the sun for a given Julian date. 


INPUTS: 

TOTJD 


JULIAN DATE 


OUTPUTS: 


0ES(6) Orbital elements of the sun; 


REFERENCES: 

SUN Is used by coverage. 








K.-.- <v-*> 
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3.5.17 SUBROUTINE TRI 

PURPOSE: 

Given either two sides and an angle of a triangle or two angles and a 
side of a triangle, an option flag, and an ambiguity flag, find the remaining 
sides and angles of the triangle. 

INPUT - OUTPUT MATRIX: 

DEFINE P. Q, R - angles (degrees) 

X, Y, Z - sides (degrees) such that 
x Is opposite P, 
y Is opposite Q, and 
z Is opposite R. 


ANGLES, SIDES (RAD) 

P, Q. Z 
P, Q, Z 
P, 0, X 
P, Y. Z 
P, Y, Z 
P, X, Y 
P, X, Y 
X, Y, Z 
X, Y, Z 
P. Q, R 
P. Q, R 

METHOD: 

Simple trigonometry. 


L 

M 

To OUTPUT: 

ANGLES, SIDES (RAD) 

R(QUAD 1 or 2), X, Y 
R(QUAD 3 or 4), X, Y 
Y(QUAD 2 or 3), R, Z 
X(QUAD 1 or 2), 0, R 
X(QUAD 3 or 4), Q, R 
Q(QUAD 1 or 4), R, Z 
0(QUAD 2 or 3), R, Z 
P(QUAD 1 or 2), Q, R 
P(QUAD 3 or 4), Q, R 
X(QUAD 1 or 2), Y, Z 
X(QUAD 3 or 4), Y, Z 


- option Flag 

- ambiguity Flag 

INPUT: 

L M 

1 1 

1 2 

2 2 

3 1 

3 2 

4 1 

4 2 

5 1 

5 2 

6 1 

6 2 
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3.5.18 SUBROUTINE VEC 


PURPOSE: 

Computes the unit vector corresponding to a given right ascension and 
declination In degrees. 

INPUTS: 

ALF — Right ascension In degrees. 

DEL — Declination In degrees. 

OUTPUTS: 

V(3) — Unit vector corresponding to ALF and DEL. 

METHOD: 

Standard trigonometry 
REFERENCES: 

No external references. 

















3.5.19 


SUBROUTINE ZERRV 


PURPOSE: 

To Initialize a real array. 

INPUTS: 

A A real array of dimension N; 

N Dimension of A; 

OUTPUTS: 

A(I) * 0.0 where 1 < I £ N; 

REFERENCES: 

This routine Is used by C0VERA6E and AAWSIM. 








3.5.20 


SUBROUTINE ZERIV 


PURPOSE: 

To Initialize an Integer array. 


INPUTS: 

IA An Integer array of dimension N; 

N Dimension of IA; 

OUTPUTS: 

IA(I) * 0 where 11 I < N; 

REFERENCES: 

This routine Is used by COVERAGE and AAWSIM. 
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actions are utl11zed In the program 


the cosine of an angle, find the angle In degrees 
an angle In degrees, output the angle such that 0 


an angle tn degrees, find the cosine of the angle 


mm the angle In radians such 

tO PI. ^ 

■ In nattflti* flwi the angle such .that 0 £ 


Of the angle 







